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Introduction 



1.1 Four fundamental forces 

Modern physics believes that matter is held together through four fundamental 
forces, gravity, the weak interaction, the electromagnetic force, and the strong 
interaction (cf. Tab. II. ip . In this section, I briefly discuss some aspects and 
differences of these forces. 

The force with the smallest relative strength is gravity, which couples two particles 
together via their mass, and which acts on all known particles. It is the most 
"common" force in our everyday life, and describes how the apple falls down 
from the tree. Newton's theory of gravity (17th century) describes this force 
as a long-range interaction between two bodies. In modern physics interactions 
are described as an exchange of virtual particles, which "carry" the force from 
one particle to another. Such a modern theory of gravity is Einstein's theory 
of general relativity (20th century). The exchange particle for the gravity is the 
graviton with vanishing mass. The range R of the force can be estimated by the 
mass m of the exchange particles [HKS95 

R ~*L' (L1 > 
where h = 6.626 0693 x 10" 34 Js is Planck's constant, and c = 299 792 458 ms" 1 
the velocity of light. The mass of the graviton is zero and therefore the range 
of gravity is infinity, as one expects. As mentioned above, gravity is very weak 
compared to the other three forces, thus it does not play an important role in 
microscopic processes. Indeed, gravity is the only one of the four forces where 
the interaction between two particles with equal charge is attractive, and not re- 
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range [m] 


exchange 


interaction between 




strength 




particle 


two particles 
with equal charge 


gravity 


10 -38 


00 


graviton 


attractive 


weak interaction 


10~ 5 


io- 18 


W ± , Z° 


repulsive 


EM interaction 


10- 2 


00 


photon 


repulsive 


strong interaction 


1 


10- 15 


gluon 


repulsive 



Table 1.1: The four fundamental forces of nature by comparison. 



pulsive. This fact, together with the infinite range of gravity, leads to its extreme 
importance in all astrophysical phenomena. 

Another fundamental force is the so called weak interaction, which acts on all 
fermions (particles with spin 1/2). The exchange particles of the weak interaction 
are the very massive Z° and bosons, their masses are ~ 80 times larger than 
the mass of the proton |EHO + 04j . Therefore this interaction is of short range 
R « 10~ 18 m, which is ~ 1000 times smaller than an atomic nucleus. It is very 
important for all decay processes (involving fermions), e.g., the beta decay of the 
neutron 

n — > p + e~ + z7 e , (1.2) 

where n is a neutron, p a proton, e~ an electron, and z7 e an anti-neutrino. 
The electromagnetic (EM) force acts on all particles with electric charge. Maxwell 
was the first (19th century) who could unify the electric and the magnetic in- 
teractions in one theory. A modern description of this interaction is quantum 
electrodynamics (QED) (20th century). The mass of the exchange particles, the 
photons, is zero, and therefore the range of the electromagnetic interaction is 
infinity. It is, in addition to gravity, very common in our everyday life, and is 
used to light our rooms and to power our hi-fi systems. In the 20th century, 
Glashow, Weinberg, and Salam worked out that the electromagnetic and the 
weak interaction can be unified in one theory, the electroweak theory. 
Last but not least, the strongest of the fundamental forces is the strong inter- 
action, all particles which experience this force are called hadrons. The modern 
theory of strong interaction is quantum chromodynamics (QCD) (this name can 
be traced back to the name of the QCD charge: colour). The exchange particles 
are 8 different massless gluons, thus one would expect, from Eq. (jl.lj) . that the 
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range of this force is infinity. However, this is not the case, the range of the strong 
interaction is R ~ 10 -15 m (of the order of an atomic nucleus). The reason for 
this is that the gluons can interact with themselfs, which leads to the finite range 
of the force and many other complications. QCD is the theory that this thesis is 
based on, and I will focus on it in the next sections. 



1.2 Quantum Chromodynamics 

In this section I discuss some aspects of QCD which are important for my thesis, 
especially the chiral and the Z(N C ) symmetry. Certainly, this introduction is far 
from complete, for more details see e.g. |Wei95at IPS95| IWei95bt Kug97[ ICL84[ 

EssnjEHsai. 

As mentioned above, modern physics believes that QCD is the best theory we have 
to describe the strong interaction between the bosonic gluons and the fermionic 
quarks. Like all fundamental forces (apart from gravity) it is based on a lo- 
cal gauge symmetry, the SU(3) C colour symmetry. As shown in tHV72 3 it is 



a renormalizable theory, which means that divergences can be regularised via 
the introduction of a finite number of counter terms to the Lagrangian of the 
theory. Another important fact is that QCD is an asymptotically free theory 
[G W731 rPo!73j . i.e., the interaction between gluons and quarks becomes weaker 
with increasing energy. One expects that this leads to a so-called quark-gluon 
plasma (QGP) at large temperatures and/or chemical potentials |CP75j . which 
means that the quarks and the gluons are no longer confined together but can 
behave as free particles in a plasma. The phase transition between hadronic mat- 
ter, which is the state of matter under "normal" conditions in our universe today, 
and the QGP is extensively discussed in this thesis. 
The QCD Lagrangian is given by 



£ = ${vfDn - m)i> - \F^F a ^ (1.3a) 



where 

D^dr-igAp*, (1.3b) 
is the covariant derivative, and 



Fr = - + 9fa bc AtA», 



(1.3c) 
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the gluonic field strength tensor. For N c colours and Nf quark flavours, ip is a 
47V c ./Vj-diniensional spinor of quark fields, if> the Dirac conjugate spinor, / y fl the 4- 
vector of Dirac matrices, m the quark mass matrix, g the strong coupling constant, 
A£ the gluonic fields, and f a bc the structure constants of the local SU(N C ) C colour 
symmetry. As mentioned above, the QCD Lagrangian foots on this symmetry, 
which means that it is invariant under transformations of the fields C(if>, A^) = 
£(^ip,^A u a tt), where 



Q(X) = exp 



i 

. a=l 



;i-4) 



describes local SU(3) C transformations. The T a 's are the generators of the sym- 
metry group, and a a are the corresponding parameters. This transformation of 
the fields is called a gauge transformation. Note that, in contrast to global trans- 
formations, cf. Eq. (|1.8|) . a local transformation depends on the space-time point 
X. In most cases I suppress the argument of Q for simplicity. The matrices Q 
satisfies 

0*0 = 1, detO = l. (1.5) 



The chiral U{Nj) r x U(Nf)i symmetry 



For massless quarks (m = in Eq Jl.Haj) the Lagrangian is invariant under the 
global chiral symmetry. To see this, the spinor and the conjugate spinor of the 
quark fields are decomposed into right- and left-handed spinors 



if) -> lf) r + lf)£ E 

Using this, the Lagrangian becomes 



1 + 75 , . 1 - 75 , 

— « — v ^ — ^ — v- 



;i.6) 



;i.7) 



Invariance under chiral symmetry means that the Lagrangian does not change, if 
the spinors are transformed via a global U{Nf) r x U(Nf)i symmetry transforma- 
tion, C(ip r , ipe) = £(Q r ip r , Qeif>e), where the right- and left-handed transformations 
are given by 



Q r = exp 







f "f- 1 \ 


• E < T a 


, = exp 


* E <AT a 


v a=0 / 




\ a=0 J 
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a r £ are the parameters, and T a the generators of the U(Nf) r ^ symmetry group. 
The chiral symmetry group is isomorphic to the vector and axial vector group, 
U(N f ) r x U(N f ) t = U(N f ) v x U(N f ) A , with V = r + £, A = r-i. Every unitary 
group can be decomposed into a direct product of a special unitary group and a 
complex phase, U(N) = SU(N) x U(l), hence the chiral group can be written as 
U(Nf) r xU(Nf)e = SU(Nf) r xSU(N f ) e xU(l) r xU(l) e ^ SU(N f ) v xSU(N f ) A x 
U(l) v x U(1) A . 

The U(l)v symmetry corresponds to baryon number conservation, and is always 
respected. In the vacuum, a non-vanishing expectation value of the quark con- 
densate, (qeq r ) ^ 0, spontaneously breaks the above symmetry to SU(Nf)v- 
This gives rise to Nj Goldstone bosons which dominate the low-energy dynam- 
ics of the theory. As shown by 't Hooft |tH76a[ ltH76bj . instantons break the 
U(1) A symmetry explicitly to Z(Nf) A |PW84aj . (For the low-energy dynamics 
of QCD, however, this discrete symmetry is irrelevant.) Consequently, one of the 
Nj Goldstone bosons becomes massive, leaving Nj — 1 Goldstone bosons. 
The SU (Nf) r x SU (Nf)gX U(1) A symmetry of the QCD Lagrangian is also explic- 
itly broken by nonzero quark masses. The Nj — 1 low-energy degrees of freedom 
then become pseudo-Goldstone bosons. For M < Nf degenerate quark flavours, 
a SU(M)v symmetry is preserved. In nature the quark masses are not equal to 
zero, i.e., the chiral symmetry is only approximatively conserved. Normally, one 
uses models which foot on the chiral symmetry only to describe the light quarks. 
The two lightest quarks are the up and down quark, 1.5 MeV < m u < 4.5 MeV 
and 4 MeV < rrid < 8 MeV. Treating these quarks as massless, i.e., using a 
SU{2) r x SU(2)i = 0(4) symmetry, yields six Noether currents, 

Jh = $r-fj1>r, Jta = $tfj1>i, (l-9a) 

with 

<Vr = dpJt = (1.9b) 

where r a are the three Pauli matrices jlZ85j . It is common to introduce the vector 
and the axial-vector currents 

V» = J» a + J£ a = $l^, A^ = J r ^-J^ a = ^5j. (1.10) 
The corresponding conserved charges are 

Q v a = J d 3 x^(x)^(x), Q A a = J d 3 x^(x) l5 ^(x). (1.11) 
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To extend this, one can consider also the mass of the strange quark as "small" , 
80 MeV < m s < 130 MeV. Then one has to replace the three Pauli matrices by 
the eight Gell-Mann matrices A a [IZ85J. The masses of the charm, bottom, and 
top quarks are very large, 1.15 GeV < m c < 1.35 GeV, 4.1 GeV < mj < 4.4 GeV, 
and m t = 174.3 ±5.1 GeV, hence introduction of these degrees of freedom into a 
chirally symmetric model is not reasonable. 

In section 1.4, I present the linear cr-model with O(N) symmetry, which foots on 
the chiral symmetry for two quark flavours. 



The Z(N C ) symmetry 

As mentioned above, QCD is based on a local SU(3) C colour symmetry. In 
the following I show how to construct from this the global Z(3) symmetry. A 
particularly simple element of the SU(N C ) C symmetry is a constant phase times 
the unit matrix, 

fi c = exp(z0)l. (1.12) 

The constraint 

det Q c = det[exp(z0)l] = exp(ziV c 0) = cos(iV c 0) + % sin(iV c 0) = 1 (1.13) 

leads to a restriction of the parameter 

2n 

<f> = Jf3> J = 0,l,...,(iV c -l). (1.14) 

Hence fl c is just the N c th root of unity times the unit matrix. Since j is an 
integer, Q c describes a global Z(N C ) rotation (symmetry) of the fields. Indeed, 
we will see that this is not a symmetry of the theory with quarks, because this 
transformation violates the periodicities of the quark fields. 
To introduce the temperature T, I work in Euclidean space-time and use the 
imaginary-time formalism |Das97[ Kap93[ ILBOOj . From the laws of quantum 



statistic, cf. |G.PY81j . we know that the (bosonic) gluon fields A u a have to be 
periodic in (0, 1/T), 

A" a (I, x) =i4£(0,x), (1.15) 



and the (fermionic) quark spinor and conjugated spinor, if> and ■0, anti-periodic, 
V>(~ x)=-V>(0,x), ffixj-^x). (1.16) 
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Obviously, any gauge transformation of these fields which are periodic in (0, 
i.e., f2(x, l/T) = fi(x, 0), does not change these periodicities, 't Hooft |tH78j was 
the first to notic that one can also consider a gauge transformation which is only 
periodic up to a constant equal to the identity matrix times the N th root of 
unity: 

n(x,^=n B , n(x,o) = i. (i.i7) 

Using this transformation the periodicity of the gluon fields remain 

nlK (x, I) n c = a: (x, I) = <(x, o) (i.i8) 

where one uses the fact that Q c , as a constant phase times a unitary matrix, 
commutes with the fields A v a , and the relation fijf2 c = 1. But the periodicities of 
the quark fields are violated 

^(x,i) = expH0)v(x,i) ^(x,0), (1.19a) 
n$(x,i) = exp(-z^(x,~) ^(x,0). (1.19b) 

Hence, the Z(N C ) is an exact symmetry only for the pure gauge theory, i.e., 
without quark fields included. In section 1.4, I present a theory which is based 
on this symmetry, the Polyakov-loop model. 



1.3 Phase transitions in QCD 

In this section I discuss the possible phases of strongly interacting matter and the 
transitions between them. Therefore, I will briefly motivate some major concepts 
for the well-known example of water. 

Thermodynamic systems can show different types of macroscopic behaviour (or in 
other words can be in different phases), depending on external degrees of freedom 
(e.g. temperature T, pressure p, chemical potential /i). Water can be in three 
different phases, the liquid phase, the gas phase, and the solid phase. By changing 
an external degree of freedom the system can change from one phase into another. 
Let's consider the following example of our everyday life. If one puts a pot with 
(liquid) water on the stove and switches the latter on, the temperature of the 
water increases up to a (critical) temperature T c = 100° C, where the water 
becomes vapour (changes the phase). The order of a phase transition is defined 
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Figure 1.1: A sketch of the phase diagram of water in the pressure vs. temperature 
plane. 



via the properties of the pressure, and its nth derivatives with respect to the 
temperature and/or chemical potential, 



d n p 



g T r. 



. d n p 
and 



dfi n 



(1.20) 

T=const. 



/x=const. 

If the pressure is continuous, but there is a discontinuity in the first derivatives, 
we call this a first-order phase transition. If the pressure and its first deriva- 
tives are continuous but there is a discontinuity in its second derivatives, we call 
this a second-order phase transition. In principle, one can define even higher- 
order phase transitions, but in nature only the first- and second- (maybe third- 
|GW801 ILM05bl ILM05al ILSSSlj ) order phase transitions occur. If there is no 
discontinuity, but a rapid change in the thermodynamic quantities, we call this a 
crossover transition. 

An appropriative way to visualise the possible phases and phase transitions of a 
system is a phase diagram. In Fig. 11.11 a sketch of the phase diagram for water 
is shown in the pressure vs. temperature plane. The lines represent first-order 
phase transitions (directly on the line both phases coincide). The first-order 
phase transition between the liquid and the gas phase becomes a second-order 
phase transition at a critical point (p, T) = (221 bar, 374 °C), and for even higher 
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pressures p > 221 bar a crossover. Another remarkable feature of this diagram is 
the triple point at (p,T) ~ (0.006 bar, 0.01 °C) where all three phases coincide. 
The transition of the example with the pot of water is a first-order transition at 
(p,T) = (1 bar, 100 °C) in this diagram. In the next sections, we will see that 
most of the discussed features can also be found in the QCD phase diagram. 

The quark mass dependence of the phase transition 

An important parameter for the order of the phase transition in QCD is given 
by the mass of the quarks [m in Eq. (jl.3a|) ]. As mentioned in Sec. 11.21 nonvan- 
ishing quark masses break the symmetry of the QCD Lagrangian explicitly. The 
principle of universality states that a phase is characterised by its macroscopic 
behaviour, which is mostly driven by the global symmetries. Therefore one way 
to find the order of the transition under this broken symmetry is to construct 
a model which respects the same (broken) global symmetry as QCD, and calcu- 
lates the order in this model instead of in QCD. In this section I consider the 
case of vanishing chemical potential \x = and three (possible) quark flavours 
Nf = 3. Therefore, the mass parameters are the strange quark mass m s and the 
degenerate masses of the up and down quark m ud = m u = m d . 
The expected phase diagram in the plane of these two mass parameters is shown 
in Fig. 11.21 This diagram exhibits two areas of first-order phase transitions, 
near m s = m u d = and m s = m u d = oo, and a crossover regime between 
them, separated by lines of second-order phase transitions. The question marks 
indicate the expected position for the physical quark masses in this diagram. If 
one of the mass parameters is infinity (i.e., the quark is infinitely heavy) the 
corresponding quark flavours are thus removed from the spectrum of physical 
excitations. Therefore one can distinguish between four different flavour cases. 
First the pure gauge theory (upper right corner of Fig. 11.2)1 where m s = m ud = oo, 
in this point all quarks are removed from the theory and only the gluons remains. 
Second, the one-flavour case Nf = 1 (right boundary of Fig. 11.2)1 where m s < oo 
but m u d = oo, on this line only the strange degree of freedom remains in the 
theory. Third, the two-flavour case Nf = 2 (upper boundary of Fig. 11.2)1 where 
m s = oo but m U( i < oo, here the strange quark is removed from the theory and 
only the up- and down quarks remain. And finally the full three-flavour case 
Nf = 3 (the rest of the plane in Fig. 11.2)1 where m s < oo and m u d < oo, in this 
area all quarks are included in the underlying theory. 
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Figure 1.2: The expected phase diagram in the plane of strange vs. degenerate up and 
down quark masses \LP03}. 

One can use universality arguments to determine the order of the phase transition. 
According to universality, the order of the chiral transition for vanishing quark 
masses in QCD is identical to that in a theory with the same chiral symmetries 
as QCD. This argument was employed by Pisarski and Wilczek jPW84a| who 
found that for Nj = 2 flavours of massless quarks (upper left corner of Fig. 11.2)1 . 
the transition can be of second-order, if the U(1)a symmetry is explicitly broken 
by instantons. It is driven first-order by fluctuations, if the U(1)a symmetry is 
restored at the critical temperature T c . For Nf = 3 massless flavours (lower left 
corner of Fig. [i~2")) . the transition is always first-order. In this case, the term which 
breaks the U(1)a symmetry explicitly is a cubic invariant, and consequently drives 
the transition first-order. In the absence of explicit U(1)a symmetry breaking, 
the transition is fluctuation-induced of first-order. 

For nonzero quark masses, the chiral symmetry of QCD is explicitly broken. 
Nonzero quark masses act like a magnetic field in spin systems, such that a 
second-order phase transition becomes a crossover transition. When the quark 
masses increase, a first-order phase transition may for a while remain of first- 
order, but it will ultimately become a crossover transition, too. In order to 
decide whether this happens for a particular choice of quark masses, universality 
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Figure 1.3: A sketch of the QCD phase diagram in the temperature vs. chemical 
potential plane. 



arguments cannot be applied, and one has to resort to numerical calculations. 
For the pure gauge theory (the theory without quarks, upper right corner of 
Fig. EU) the transition is of first-order [B+881 IBCD+88I IUka901 IUka971 ILae981 
IKarOOj as predicted by Svetitsky and Yaffe [SY82 . Nevertheless, recent results 
show that the transition is a weak first-order transition, thus it is more accurate 
to speak of a "nearly second-order" transition (KKLLOOj . 

The temperature vs. chemical potential plane 

In Fig. 11.31 a schematic plot of the QCD phase diagram in the temperature vs. 
chemical potential plane is shown. This diagram exhibits mainly three different 
phases, the hadronic phase, the QGP, and the colour-superconductivity phase. In 
the hadronic phase, the energy is small enough for stable hadrons to exist, e.g., 
protons, neutrons, pions, etc. As mentioned above, QCD is an asymptotically free 
theory, therefore for large temperatures and/or chemical potential the quarks in 
the hadrons become deconfined and behave more or less as free particles in a 
plasma, the QGP. There is probably a third phase at small temperatures and 
very large chemical potential, the phase of colour superconductivity, where a 
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quark-quark Cooper pair can be build, in analogy to the well-known electron- 
electron Cooper pair of electromagnetic superconductivity. Note that this phase 
can be separated into many other phases, e.g., the 2SC, CFL, CSL, or the polar 
phase, c.f. |BL84l IA RW98I IRSS V981 [AlfflTI iRMMl ISho04l iRWOOj . 
First, let's consider the transition between the hadronic phase and the QGP. 
The line between these two phases corresponds to a first-order phase transition 
which ends in the chiral critical endpoint of a second-order phase transition. The 
position of this chiral critical endpoint is still under investigation l'K02i [FK04, 
ISte04l IGGG05I IBCPC94I [AY891 ITTT031 [A1T031 lE+04] . the latest lattice results 
from Fodor and Katz [FK04 found this point at T = 162 ± 2 MeV and fi = 
120 ± 13 MeV (indeed at nonzero chemical potential lattice QCD calculations are 
hampered by the fermionic sign problem and become increasingly unreliable with 
increasing chemical potential |dFP02j ). As indicated in this plot, there are two 
possibilities for an order parameter of this phase transition, the quark condensate 
(qq) which links the transition to the restoration of the chiral symmetry, and 
the expectation value of the Polyakov-loop (/) which links the transition to the 
restoration of the Z(N C ) symmetry. In contrast to the quark condensate, the 
expectation value of the Polyakov-loop is zero in the hadronic phase and nonzero 
in the QGP, i.e., the Z(N C ) symmetry is broken in the QGP and restored in 
the hadronic phase. However, note that these order parameters correspond to 
two very different Ansatze for the underlying symmetry. The Z(N C ) symmetry 
is exact for infinitely heavy quarks (in the pure gauge theory), and the chiral 
U(Nf) r x U(Nf)i symmetry for quarks with vanishing mass. Although there 
are lattice QCD calculations which indicate that these two Ansatze lead to the 
same phase transition !Kar02j , the relationship of the symmetries is still an open 
question in modern high energy physics, cf. e.g. jMST04| ISKT04| IKL99 . 
The stable ground state of nuclear mater at vanishing temperature is at nonzero 
chemical potential fiQ = 308 MeV. For smaller chemical potential, nuclear matter 
shows a similar behaviour to a gaseous phase and for larger the behaviour of a 
liquid. As the line of first-order phase transitions in the phase diagram of water, 
this line in the QCD phase diagram ends in an critical point of second-order (at 
a temperature of ~ 10 MeV). 

Note that the area of colour-superconductivity is still not well understood, espe- 
cially for "small" chemical potential. Maybe there is no direct phase transition 
between the hadronic phase and the colour superconductivity phase without pass- 
ing the QGP, i.e., maybe there is no triple point between the three phases (indi- 
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cated by the question marks in Fig. 1.3). Or maybe there is no phase transition 
at all |ABR99[lRW00j . 

1.4 Effective models of QCD 

At temperatures of the order of ~ (qq) 1 ^ 3 , the thermal excitation energy is large 
enough to expect the restoration of chiral symmetry. At such energy scales, 
the QCD coupling constant is still large, rendering perturbative calculations un- 
reliable. Thus, one has to resort to nonperturbative methods to study chiral 
symmetry restoration. A first-principle approach is lattice QCD |Kar02j . Lat- 
tice QCD calculations have determined the temperature T c for chiral symmetry 
restoration to be of order 150 MeV at zero quark chemical potential [LP03 : . These 
calculations, however, face several technical problems. 

The first is that they become numerically very difficult for physically realistic, 
i.e., small, values of the up- and down-quark masses. Although progress in this 
direction has been made jFK04| . most studies use unphysically large values. An- 
other problem is that, at nonzero chemical potential, lattice QCD calculations 
are hampered by the fermion sign problem and become increasingly unreliable 
for chemical potentials larger than (a factor tt times) the temperature }dFP02j . 
An alternative nonperturbative approach to study chiral symmetry restoration 
is via chiral effective theories. These theories have the same global U(Nf) r x 
U(Nf)e symmetry as QCD but, since quark and gluons are integrated out, do 
not possess the local SU(3) C colour symmetry of QCD. The effective low-energy 
degrees of freedom are the (pseudo-) Goldstone bosons of the QCD vacuum, i.e., 
the pseudoscalar mesons. However, in the chirally symmetric phase these particles 
become degenerate with their chiral partners, the scalar mesons. Therefore, an 
appropriate effective theory to study chiral symmetry restoration in QCD is the 
linear sigma model |Lev67t IGML60| which treats both scalar and pseudoscalar 
degrees of freedom on the same footing. The advantage of chiral effective theories 
over lattice QCD calculations is that their numerical treatment (within some 
many-body approximation scheme) is comparatively simple and that there is no 
problem to consider arbitrary quark chemical potential. 

In the following, I introduce a reasonable effective model which respects the chiral 
symmetry, the linear a- model with O(N) symmetry, and an effective model which 
is based on the Z(N C ) symmetry, the Polyakov-loop model. 
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The linear cr-model with 0(N) symmetry 
The Lagrangian of the O(N) linear sigma model is given by 

C(4>) = ld lt <f>-&'4>-u(4>) (1.21) 

where cf> = (0i,7r), with the first component 0! corresponding to the scalar 
cr-meson and the other components tz = (0 2 , ...,4>n) corresponding to the pseu- 
doscalar pions. (However, note that the original cr-model introduced by Gell- 
Mann and Levy [GML60 , IZ85 incorporates the scalar and pseudoscalar mesons 
as well as a fermionic isodoublet field if) of mass zero. In the linear cr model 
discussed in the following, the fermionic degrees of freedom are neglected.) The 
function U(<f>) is the tree-level potential, 

U{4>) = \l?<t> ■ + ^(0 » 2 - Hfa , (1.22) 

where /i 2 is the bare mass and A > the four-point coupling constant. For 
fi 2 < the O(N) symmetry is spontaneously broken to 0(N — 1), leading to 
N — 1 Goldstone bosons, the pions. The parameter H breaks the symmetry 
explicitly, giving a mass to the pion. 

The parameters can be expressed in terms of the vacuum mass of the cr-meson, 
rricr, the vacuum mass of the pion, m % , and the vacuum decay constant of the 
pion, f v , 

2 _ ml -3ml _ N(ml - m\) 2 f , 

V ~ 2 ' 8^ ' H ~ m ^^- t 1 - 26 ) 

With m a = 600 MeV, m w = 139.5 MeV, and / ff = 92.4 MeV jEHO+04j this leads 
to H = (121.60 MeV) 3 , A = 19.943, and /i 2 = -(388.34 MeV) 2 . 
I assume translational invariance so that one may consider the effective potential 
V instead of the effective action T, cf. Sec. 1.5. For translationally invariant 
systems, these two quantities are related via 

r[cr, 7T, S, P] = V[a, 7T, S, P] , (1.24) 

where Sl 3 is the 3-volume of the system, and cr, it, S, P are the expectation values 
of the one- and two-point functions for the scalar and pseudoscalar fields in the 
presence of external sources |G,IT74j . I am interested in the case where these 
sources are zero, i.e., in the stationary points of F or V. Because the vacuum of 
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Figure 1.4: The set of two-particle irreducible diagrams considered in the linear a model 
with O(N) symmetry. The diagrams a,b, and c are the double-bubble diagrams, and d 
and e are the sunset diagrams. A full line denotes the full propagator for the a-meson 
and a dashed line the full propagator for the pion. The four-particle vertex ~ A is 
represented by a square and the three-particle vertex ~ Act by a triangle. 



QCD has even parity, the expectation values of the pseudoscalar fields are zero, 
fr = 0, and I shall simply omit the dependence of V on fr in the following. 
Then, the effective potential for the O(N) model in the CJT formalism, discussed 
in Sec. 1.5, reads |RRR031 ILROO] 



V[a,S,P] = U(a) + - \nS- 1 (Q) + S- 1 (Q;a)S(Q)-l 

2 Jq l 



+ 



N-l 



\nP- 1 (Q)+p- 1 (Q;a)P(Q) - 1 



+ V 2 [a,S,P] 



;i.25) 



where U(a) is the tree-level potential ()1.22|) . evaluated at 4> = (a, 0, . . . , 0). At 
tree-level the one-point expectation value can be calculated analytically 



a 



Nli 2 2 e „ 

= cos -, a — arccos 

4A v/3 3' 



HN 
~8A~ 



12A 



3/2' 



(1.26) 



The quantities S 1 and P 1 are the inverse tree-level propagators for scalar and 
pseudoscalar mesons, 



S~\K-a) 
F-\K-a) 



-K 2 + ml(a) 



(1.27a) 
(1.27b) 
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where the tree-level masses read 



ml(a) 



o 12 A o 

2 -2 



(1.28a) 
1.28b) 



As explained in Sec. 11.51 the functional V2 in Eq. ()1.25|) is the sum of all 2PI 
diagrams. The standard Hartree-Fock approximation is defined by restricting 
this sum to the three double-bubble diagrams shown in Figs. 11.41 a. b, and c, 



v 2 a [P] = (iv + i)(iv-i)A 



HQ) 



vm = 3 a 



S(Q) 
A 



V 2 C [S,P] = 2(N-l)-[ S(Q)J l P(L) 



(1.29a) 
1.29b) 
(1.29c) 



Later, we will see that, in order to include the nonzero decay width of the par- 
ticles, one has to go beyond the Hartree-Fock approximation by additionally in- 
cluding the sunset diagrams of Figs. I1.4l d and e. These diagrams have an explicit 
dependence on a, 



V 2 d [a,S,P] 



V£[cj,S} 



§{L)P{Q)P(L + Q) , (1.30a) 



1 3! 
2 



'4Act\ 2 



N 



S(L)S(Q)S(L + Q) . 



1.30b) 



The complete expression for V2 is obtained by the sum of the contributions in 
Eqs. (fT2HD and ([Oljl . 



v 2 = v 2 a + v 2 b + V 2 C + v 2 d + v 2 e . 



(1.31) 



The corresponding Dyson-Schwinger equations for the dressed propagators and 
the condensate equations can be derived by the minimisation of this effective 
potential, which will be the topic of the following chapters. 



The Polyakov-loop model 

In this section, I discuss how to construct an effective model, based on the Z(N C ) 
symmetry, as discussed in Section 1.2. To find an order parameter for the pure 
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gauge theory with N c colours let's first consider the Wilson line |Pol78[ |£H78, 
I7S82] 



Pexp 



ig 



l/T 



Aq(x, r)dr 



1.32) 



where A denotes the temporal component of the gauge field in the fundamental 
representation, g the gauge coupling, and V denotes path ordering with respect 
to imaginary time r. The thermal Wilson line transforms under local SU(N C ) 
transformations as: 



Lfx) 



fit ( X) i)L(x)nt(x,0). 



1.33) 



An order parameter for the Z(N C ) symmetry is given by the trace over the Wilson 
line 



N, 



tr V exp 



l/T 



A (x, t) dr 



1.34) 



Its expectation value, io(T), vanishes when T < T c , and is nonzero above T c . 
Indeed, by asymptotic freedom, £ — > 1 as T — > oo. The simplest guess for a 
potential for the Polyakov-loop is: 



V(£) 



b i\t\ 2 + -A\t\ 2 )\ (N c = 2) 



1.35) 



2 1 1 4 

The Polyakov-loop model jFjOQl lPis02l IDP011 lDPT)2bl IDP02al ISD.lOlj is a mean- 
field theory for t. In a mean-field analysis all coupling constants are taken as 
constant with temperature, except for the mass term, ~ &2^| 2 - About the transi- 
tion, condensation of £ is driven by changing the sign of the two-point coupling: 
b 2 > above T c ^(T) — > 1 for T — > oo], and < below T c . 
For two colours, p.35|) is a mean-field theory for a second-order deconfining tran- 
sition |ES99l lEWl IEMSZ961 IEKR951 lEFR+89] . The I field is real, and so the 
potential defines a mass: (m^/T) 2 = (l/Z s )d 2 V/d£ 2 , where Z s is the wave func- 
tion renormalisation constant for I |Wir02[ ID()03j . The mass is measured from 
the two-point function of Polyakov-loops in coordinate space, oc (1/r) exp(— mi r) 
as r — > oo. 

For three colours, £ is a complex valued field, and a term cubic in I appears in 
V(l), 



w 2 1 1 3 2 



+ l(l^| 2 ) ; 



(N c = 3) 



(1.36) 



At very high temperature, the favoured vacuum is perturbative, with £q w 1, 
times Z(3) rotations. One then chooses 63 > so that in the Z(3) model, there is 
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always one vacuum with a real, positive expectation value for £ . This produces 
a first-order deconfining transition, where £q jumps from at T 6 7 to t c = 2& 3 /3 
at T+ jDPnil IDPn2bl IDPn2aj : T c is given by b 2 {T c ) = -2b 2 3 /9. The £ field has 
two masses, from its real {mi) and imaginary (m^) parts. At T c + , yJ~Z~ s mi/T = £ c . 
The mass for the imaginary part of £ is \f^ s rhi(T) /T oc \fb^£] at T c + , rhi/mi = 3, 
twice the value expected from a perturbative analysis of the loop-loop correlation 
function, obtained by expanding £ from Eq. ()1.34|) to order Aq [DP 02cI1 IDP02cj . 
This mass ratio receives corrections if five-point and six-point couplings are in- 
cluded in the effective Lagrangian [DP02d, ll)P02cj but those are not crucial 
for the following discussion. One notes that, in principle, all of the above cou- 
pling constants could be determined on the lattice. The lattice regularization 
requires non-perturbative renormalization of the Polyakov-loop in order to define 
the proper continuum limit of £ |K K PZ02I iPPfMl IZan03l IDHL+04| . 
Within the above mean-field theory, dynamical quarks act like a "background 
magnetic field" which breaks the Z(3) symmetry explicitly, and a term linear in 
£ also appears in V(£) BU831 IHK8! IMQ951 11+99] : 

m = -6^ - h W _ f^f + I {W f, {Nc = 3 , m , < co) . 

(1.37) 

Hence, as m n decreases from infinity, bi(m n ) turns on. The normalisation of 

for T — > 00 is such that £q 1, i.e., 6 2 (^ = 00) = 1 — b\ — 63. 

In Chapter II, the Polyakov-loop model is used to study the QCD phase transition. 



1.5 The Cornwall- Jackiw- To mboulis formalism 

At nonzero temperature T, ordinary perturbation theory in terms of the cou- 
pling constant g breaks down in the sense that one can no longer order different 
contributions according to powers of g |D J74j . This is because the new energy 
scale introduced by the temperature can conspire with the typical momentum 
scale p of a process so that gT/p is no longer of order g, but can be of order 1 
BP90a] IBP9 0bj. Consequently, all terms of order gT/p have to be taken into 
account which requires the resummation of certain classes of diagrams. 
A convenient resummation method is provided by the extension of the Cornwall- 
Jackiw-Tomboulis (CJT) formalism |CJT74j to nonzero temperatures and chem- 
ical potentials. The CJT formalism is equivalent to the ^-functional approach 
of Luttinger and Ward |LW60j and Baym |Bay62| . It generalises the concept 
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of the effective action r[0] for the expectation value <f) of the one-point func- 
tion in the presence of external sources to that for the effective action T[4>, G] 
for one and two-point functions, and G, in the presence of external sources. 
(For an extension of this approach to three- and more-point functions, see 

The starting point for the CJT formalism is the generating functional for Green's 
functions |CJT74j 

Z[J,K] = e w[J ' K] = |p0exp|/[0] + <j, J +-</,!( </>} , (1.38) 

where J[0] = J x £ is the classical action, W[J, K] the generating functional for 
connected Green's functions, and with the short-hand notations 



(p(X)J(X), (1.39a) 
/ 4>(X)K(X,Y)<t>(Y). (1.39b) 

JX,Y 



IX,Y 

The expectation values for the one- and two-point functions in the presence of 
the external sources, and G, are given by 

**> - ^ 

Note that the derivative with respect to K is 

M(X-Y) = \ [0{X ' Y) + ^ X) ^ Y)] - (L41) 

The aim is now to eliminate the sources in favour of the one- and two point 
functions with a double Legendre transformation 

G] = W[J, K}-0J- hj>K4> - X -GK, (1.42) 

where GK = J XY G(X, Y)K(X, Y). Minimisation of this functional gives the ex- 
pectation values of the one- and two-point function in the absence of the external 
sources, <p and Q, 

srfoG] 



5<f>{X) 
5G{X, Y) 



= 0, (1.43a) 

=<p,g=g 

= 0. (1.43b) 

=v,G=g 



28 



Introduction 



The last equation corresponds to the Dyson-Schwinger equation for the full 
dressed propagator. In the CJT formalism, one uses the following effective action 
C.)T74| 



where, 



ilr HG- 1 ) 



G-\X,Y; 



s 2 m 



1.44) 



1.45) 



50(X)50(F) 

is the tree- level propagator, and r 2 [0, G] the sum of all two-particle irreducible 
(2PI) vacuum diagrams with internal lines given by G. I'm interested in the 
case where the fields are translationally invariant, so that one may consider the 
effective potential V instead of the effective action. For translationally invariant 
systems, these two quantities are related via T[4>, G] = —Qz/TVl&G] , where Q3 
is the 3-volume of the system, thus 

V[4,G\=U@) + l [inG-'ify + l [[G- 1 (k;4>)G(k)-l] + V 2 [^G], (1.46) 
where 



2 Jk 



2 Jk 



VS,G] = -TT 2 &G}/n 3 . (1.47) 

The stationary conditions in Eqs. (|1.43j) . can now be expressed in terms of the 
effective potential 

SV ^ = 0, (1.48a) 

= 0. (1.48b) 



S(f> 
5V$,G] 



5G(K) 



=V,G=Q 



Together with Eq. (jl.46|) . the last stationary condition gives the Dyson-Schwinger 
equation 

g- 1 (k)=G- 1 (k;<p) + Tl(k), (1.49) 

where 

SV 2 [0,G] 



U(k) = 2 



5G{k) 



1.50) 



is the self-energy. Note that the thermodynamic pressure is (up to a sign) identical 
to the effective potential |Riv88j 



P=-V[<p,Q\. 



1.51) 



1.6 The aim of this work 
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1.6 The aim of this work 

The main focus of my thesis is QCD and its phase transitions as discussed in the 
introduction. 

In chapter II, I present a study I've done with Adrian Dumitru and Jorg Ruppert 
DR R04j about the phase transition temperature of QCD, and its dependence 
on the quark (or pion) mass. In the first part of this chapter, I use the linear 
u-model with O(N) symmetry (cf. Sec. 1.4), within the CJT formalism (cf. Sec. 
1.5), which links the transition to chiral symmetry restoration. The parameters 
of this model depend on the vacuum mass of the pion and cr-meson and the 
vacuum value of the decay constant of the pion. From lattice QCD calculations 
we know the dependence of these values on the quark mass, and hence one can 
determine the dependence of the chiral phase transition (especially the transition 
temperature) on the quark mass. In the second part of this chapter, I use the 
Polyakov-loop model (cf. Sec. 1.4), which links the transition to the restoration 
of Z(N C ) symmetry. The aim of this part is to find how "strong" one has to break 
the Z(3) symmetry in order to reproduce the transition temperature (given by 
lattice QCD calculations). 

In chapter III, I present a study I've done with my Ph.D. advisor Dirk Rischke 
and Jorg Ruppert RRR05J about the improvement of the standard Hartree- 
Fock approximation by including nonzero decay width effects. In the standard 
Hartree-Fock approximation one treats all particles as stable quasiparticles, which 
means that the spectral densities of these particles are just delta-functions with 
zero decay width, T = 0. Obviously, this is a reasonable approximation for all 
particles with small decay width, i.e., the pion with a decay with of ~ 8 
eV in vacuum, but not for broad particles, like the a-meson with a vacuum 
decay width of T a ~ (600 — 1000) MeV. The decay width is proportional to 
the imaginary part of the self-energy of the particle T ~ ImlL In the standard 
Hartree-Fock approximation only real-valued tadpole diagrams (cf. Fig. 3.1 b, c, 
and 3.2 a, c) are taken into account. I improve this scheme by taking additionally 
into account the cut sunset diagrams (cf. Fig. 3.1 d, e, and 3.2 d), which 
have a real and an imaginary part. Beside the imaginary part, another major 
difference between these two types of diagrams is that the cut sunset diagram 
is a functional of an external four momentum vector K = (&>,k) (where to is 
the energy and k the momentum of the particle). Hence, the Dyson-Schwinger 
equations for the full propagators and the condensate equations become integral 
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equations one has to solve on an energy-momentum grid. In this chapter, my 
focus is to study the influence of the inclusion of nonzero decay width effects, 
therefore, as a first approximation, I simply neglect the real part arising from 
these cut sunset diagrams. 

Finally, in chapter IV, I present a work |Rod 05 I've done to study the effects 
of the (in chapter III neglected) 4-momentum dependent real parts of the cut 
sunset diagrams in the linear cr-model with O(N) symmetry. I do this not in 
the full Hartree-Fock but in the Hartree approximation. The difference between 
these both approximation is that all terms of order ~ 1/N are neglected on the 
level of the condensate and Dyson- Schwinger equations in the Hartree approxi- 
mation. This leads to a vanishing imaginary part of the pion self-energy, hence to 
a vanishing pion decay width. The decay width of the cx-meson remains nonzero. 
Besides the effects of the real part of the cut sunset diagram, I study in this 
chapter the influence of the "choice" of the vacuum mass of the cx-meson. As 
mentioned above, the parameters depend on the vacuum mass of the cx-meson 
mass, which is not well defined, m a = (400 — 1200) MeV, because of the very 
large (vacuum) decay width of the er-meson. Therefore, I compare the results for 
m a = 400, 600, and 800 MeV, in the case of explicit chiral symmetry breaking 
ma 7^ and the chiral limit = 0. 



1.6 The aim of this work 
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Conventions 

I denote 4- vectors by capital letters, X = (xo,x), with x being a 3- vector of 
modulus x = |x|. The imaginary-time formalism is used to compute quantities 
at nonzero temperature. Integrals over 4- momentum K — (u,k) are denoted as 

/ K /W^rJj|l/(-^,k), (1.52) 

where T is the temperature and u n = 2irnT, n = 0, ±1, ±2, . . . are the bosonic 
Matsubara frequencies. My units are H — c — ks — 1. The metric tensor is 

<T = diag(+, -,-,-). 
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The quark mass dependence of 
the transition temperature 

2.1 Motivation 

Lattice QCD calculations at finite temperature and with dynamical fermions are 
presently performed for quark masses exceeding their physical values; for a recent 
review see jTTP03] . To date, pion masses as low as ~ 400 MeV are feasible IKLPOlj , 
about three times the physical pion mass. When comparing effective theories to 
first-principles numerical data obtained on the lattice it is therefore important 
to fix the parameters (coupling constants, vacuum expectation values and so on) 
such as to match the values of physical observables, e.g. of m-^, to those of the 
lattice calculations. For example, the QCD equation of state in the confined phase 
appears to be described reasonably well by that of a hadron resonance gas model, 
after extrapolating the physical spectrum of hadrons and resonances to that from 
the lattice |KRT03I IKLPOlj . Thus, lattice dat a on the dependence of various 
observables on the quark (or pion) mass constrain effective theories for the QCD 
phase transition at finite temperature and could provide relevant information on 
the driving degrees of freedom. 

In the following, I analyse the dependence of the chiral symmetry restoration 
temperature on the vacuum mass of the pion using a linear sigma model in section 
2.2. The linear sigma model provides an effective Lagrangian approach to low- 
energy QCD near the chiral limit [PW84b, RW93J. It incorporates the global 
flavour symmetry, assuming that "colour" can be integrated out. For example, 
it allows one to discuss the order of the Nf = 2 + 1 chiral phase transition as 
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a function of the quark masses jPW84bl IH.W93I [LRSBOQl IMMOP94I IM()PP92[ 
IHoTS3ll(^P94| . 



Instead of working up from zero quark mass, one could start with the quark 
masses taken to infinity, that is, with a pure gauge theory. Then, one can discuss 
the deconfinement transition at finite temperature within an effective Lagrangian 



for the Polyakov loop with global Z{N C ) symmetry |YS82l lPis02l IDP011 IDP02H 
IDPf)2al LSD.) Oil IOM00I IMMO021 IMST04I IFuk04| (JV C is the number of colours). 



For finite pion mass, the symmetry is broken explicitly, and the phase transition 
(or crossover) temperature is shifted, relative to the pure gauge theory where 
pions are infinitely heavy. In section 2.3., I determine the endpoint of the line of 
first-order transitions for three colours, and extract the magnitude of the explicit 
Z(3) breaking from lattice data on AT C . 

Pion Mass and Decay Constant in Vacuum 

As discussed in Sec. I1.2[ the Lagrangian of QCD with the quark mass matrices 
set to zero is invariant under independent rotations of the Nf right- and left- 
handed quark fields. It exhibits a global SU(Nf) r x SU (Nf)e symmetry, leading to 
2(Nj — 1) conserved currents. Those are Nj — 1 vector currents, Vf = ■07 M Aj , /2, 
and Nj — 1 axial currents, Af = ip'y fl 'y 5 \ i ip /2, with Aj the generators of SU(Nf), 
normalised according to tr AjAj = 25^-. The SU(Nf) v subgroup of vector trans- 
formations is preserved in the vacuum |VW84| . while the SU(Nf)A is broken 
spontaneously by a non-vanishing chiral condensate (qaqL) ^ 0, leading to non- 
conservation of the axial currents. 

In reality, of course, even SU(Nf)y is broken explicitly by the non-vanishing 
quark mass matrix. Nevertheless, since m u ~ rrid, the SU(2)v symmetry is 
almost exact in QCD. The small explicit breaking of SU(2)a is responsible for 
the non-vanishing pion mass, as given by the Gell-Mann-Oakes-Renner relation 



I neglect isospin breaking effects here, and so assume that m u = m d = m q . (qq) 
denotes the sum of the vacuum expectation values of the operators urUl and 
djidi, and their complex conjugates. The proportionality constant /„- is the pion 
decay constant. It should be noted that (|2.1|) is only valid at tree level, and that 
loop effects induce an implicit dependence of both f w and (qq) on m q . For small 




(qq) 



(2.1) 



2.2 Results: The linear ct-model with O(N) symmetry 
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m q , this dependence can be computed in chiral perturbation theory [GL 89J. For 
example, at next-to-leading order, 



where M and F are the couplings of the effective theory (equivalent to m q and 
(qq)), and A3 and A4 are two renormalization-group invariant scales. These re- 
lations link the behaviour of f n to that of m n , the mass of a physical state. (In 
what follows, I use to vary the strength of explicit symmetry breaking rather 
than using directly the scale-dependent quark masses). 

More accurate results than Eqs. ()2.2al2.2b|) can perhaps be obtained by com- 
puting quark propagators for various quark masses on the lattice. Ref. [CH03J 
analysed the propagators for gauge-field configurations generated with the stan- 
dard Wilson gauge action ("quenched QCD"), using overlap fermions with exact 
chiral symmetry. They obtained a parametrisation of both and f n in terms of 
the mass m q of u and d quarks (see section 2 in |CH03j ) which allows to express 



f n as a function of m^. Their data covers an interval of 0.4 GeV < m n < 1 GeV, 
and 0.15 GeV < y/2f w < 0.22 GeV. 

2.2 Results: The linear cr-model with O(N) sym- 
metry 



In this section, I discuss chiral symmetry restoration at nonzero temperature, 
and in particular the dependence of the symmetry restoration temperature on 
the pion mass. For simplicity, I restrict myselfe here to the two-flavour case. My 
emphasis is not on the order of the transition as the strange quark mass is varied 
but rather on how the temperature at which the transition occurs (be it either 
a true phase transition or just a crossover) depends on the pion mass. Such 
dependence arises from two effects. First, of course, due to explicit symmetry 
breaking occurring when > 0. Second, due to the "indirect" dependence of 
spontaneous symmetry breaking, i.e., of the condensate (qq) resp. f n , on the pion 
mass (through pion loops, see previous section). 

In the following I use the counter-term renormalization scheme as discussed in 
section V.B in |LR00j . In this scheme a mass renormalization scale \i is introduced 




(2.2b) 



(2.2a) 
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and the couplings then depend on that scale. However, choosing 
2 _ [ m^(ln ml -I)- m%(ln ml - 1) ] 

Hren - eX P 2 _ 2 ' 

<T 7T 

the four-point coupling A(/i) = A tree retains its tree-level (classical) value [LROO . 
In other words, this renormalization prescription evolves the renormalization scale 
H ren in such a way as to keep A constant. 

Explicitly, this leads to the following expressions for the couplings {LROOJ: 

1 2 2 
1 m a — 

A 



ml 



P 

J 7T 



H = U {ml - 2Xf; 



where 



(4tt) 2 



M 2 In 



m 



M 2 



-\ (ml - 3m 2 n ) - 6XQ 



M 2 + /i 



2 

ren 



(2.4) 



(2.5) 



These equations determine the couplings in vacuum in terms of m n , f n and m a . 
The dependence of /„- and m n on the quark mass is taken from the data of 
ref. [CH03 . Roughly, for m w : 0.4 GeV — > 1 GeV, f n increases by about 50 %, 
leading to an increase of the explicit symmetry breaking term H by a factor of 10. 
I also require the dependence of m a on m n , which I take from a computation with 
standard Wilson fermions |K + 04bj . Those authors find that m a is essentially a 
linear function of m\. I checked how my results in Fig. 12.11 depend on this 
assumption by using, alternatively, a linear dependence m a = m n + const., with 
m a = 0.6 GeV for m, = 0.14 GeV. I found essentially the same dependence of T c 
on m n . 

At nonzero temperature, one uses the effective potential for composite operators, 
as discussed in Sec. 11.51 to determine the masses and the scalar condensate in the 
Hartree-Fock approximation. This approximation is defined by only taking into 
account the double-bubble diagrams shown in Figs. 11.41 d and e in the effective 
potential V. The gap equations for the condensate and the masses are given by 
minimisation of this potential. The expectation values of the one- and two-point 
functions in the absence of external sources, a and S, V, are determined from 
the stationary points of V, 
5V 

Id 



6V 
°<6S 



sv 



0. 



(2.6) 

a=cr,S=S,P=V WKJ a=cr,S=S,P=V WJ - 5=cr,S=S,P=T 

leading to the gap equations for the chiral condensate a, and the effective masses 
of the (j- me son and the pions, M a and M n , 



H 
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~Ml - 2Xa 2 



(2.7a) 
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Figure 2.1: Left: The crossover temperature T c as a function of the (vacuum) pion 
mass as obtained from the linear sigma model with 0(4) symmetry in comparison to 
lattice data \KLP01\ for two and three flavours. The scale for both T c and m n is set by 
the zero-temperature string tension in the pure gauge theory, ^/cx ~ 0.425 GeV. Right: 
The scalar condensate, {a), as a function of temperature for various pion masses. 



Ml = m 2 + 3A [a 2 + \Q T {M a ) + Q^M a )} + [Q T (M^) + Q M (M^)]} , (2.7b) 
Ml = m 2 + \{a 2 + [Q T (M a ) + Q ll (M a )} + 5[Q T (M 7T ) + Q fl (M 7r )}}, (2.7c) 

where the nonzero-temperature contribution of the tadpole diagram is given by 

1 r°° 

Q T (m) = — dqiouiq)}- 1 f[cu(q)}, (2.8) 

Z7T Z JO 

where f(cu) = l/[exp(o;/T) — 1] is the Bose-Einstein distribution function and 
ui(q) = y/q 2 + m' 2 the quasiparticle energy. The self-consistent solution of the 
above gap equations for a given vacuum pion mass determines the temperature 
dependence of the scalar condensate as the order parameter of chiral symmetry 
restoration. For explicitly broken chiral symmetry, H > 0, the transition in this 
approach is a crossover. I define the crossover temperature T c by the peak of 
da/dT. The dependence of T c on m n is depicted in Fig. 12.11 (left), where I have 
also shown lattice results obtained with two and three degenerate quark flavours, 
respectively jKLPOlj (the Nf = 2 data with standard action, the Nf = 3 data 
with improved p4-action). Driven by the increase of both f„ and H with m n , the 
linear sigma model predicts a rather rapid rise of T c with the pion mass, as com- 
pared to the data which is nearly flat on the scale of the figure. While lattice data 
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indicate a rather weak dependence of T c on the quark mass (see also ref. |B+97j). 
models with spontaneous symmetry breaking in the vacuum naturally predict a 
rather steep rise of T c with the VEV <r vac = /„., which itself increases with the 
quark (or pion) mass. My findings here are in qualitative agreement with those 
from ref. |BJ W99j who employed nonperturbative flow equations to compute the 
effective potential for two-flavour QCD within the linear sigma model. They also 
find a steeper slope of T c (m n ) than indicated by the lattice, even though their 
analysis appears to predict a somewhat weaker increase of f n with m n than the 
data of [CH03J, which I employ here. 

Fig. 12 .11 also shows the temperature dependence of the a condensate (right). With 
m a a linear function of |K + 04bj . the width of the crossover is approximately 
independent of the pion mass for 0.4 GeV < m n < 1 GeV, while I found consid- 
erable broadening when m a is linear in (not shown). The chiral susceptibility 
da/dT at its maximum is ~ 0.25, i.e., the crossover is in fact quite broad for the 
range of m n considered. Since this is at variance with lattice data on QCD ther- 
modynamics (pressure and energy density as functions of temperature, see e.g. the 
review in |LP03j ) , one might argue that the crossover is in fact not driven by the 
order parameter field but by heavier degrees of freedom .K RT03| IKLP01[ RTL89 . 
Such degrees of freedom could reduce the pion-mass dependence of the tran- 
sition substantially: using three-loop chiral perturbation theory (i.e., the non- 
linear model), Gerber and Leutwyler find jGL89l) that T c increases rapidly from 
« 190 MeV in the chiral limit (using their set of couplings) to ~ 240 MeV for 
physical pion mass. However, when heavy states are included (in the dilute gas 
approximation), then T c increases less rapidly, from ^170 MeV in the chiral limit 
to ~ 190 MeV for physical pion mass. 
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2.3 Results: The Polyakov-loop model 

In this section I consider the Polyakov-loop model as introduced in Sec. 1.4. As 
discussed, the effective potential for three colours with explicit breaking of the 
Z(3) symmetry is given by, 



m = - h l -±f - | w . - + 1 0*')" • - 3 - - < ») • 

(2.9) 

First, I consider the case where bi is very small, and take the term linear in I 
as a perturbation; then the weakly first-order phase transition of the pure gauge 
theory persists (in what follows, the critical temperature in the pure gauge theory 
with b\ = will be denoted T*). The critical temperature is determined from 

62(r e ) = ~6g(l + y|)+0(6?). (2.10) 
The order parameter jumps at T c , from 



Zo(T-) = l b ^ + 0(bl), (2.11) 



to 

UTr) = \h-^ + 0(b\) . (2.12) 

Note that numerically £q(T~) could be much larger than b\ if the phase transition 
in the pure gauge theory is weak and so the correlation length £ = 1/me near 
T c is large (i.e., if 6 3 is small), as indeed appears to be the case for N c = 3 
colours [KKLL00 . In other words, it could be that on the lattice I quickly 
develops a non-vanishing expectation value at T~ already for rather large quark 
(or pion) masses, but this does not automatically imply a large explicit symmetry 
breaking (see also Fig. 12.21 below). 

From Eq. ()2.10j) one can estimate the shift of T c induced by letting < oo. 
Writing the argument of b 2 in that equation as T* + AT C and expanding to first- 
order in AT C one obtains 




~UT-)hT^\ +0(bl). (2.13) 



3 \ dT 



T=T* 
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The shift in T c with decreasing pion mass is proportional to the expectation value 
of the Polyakov-loop just below T c ; all other factors on the right-hand side of 
Eq. (j2.13J) do not depend on 61 or m w . Numerical values for 6 3 and for b 2 (T) were 
obtained in |l)P01llHF02HlSl).IMllBnU12] by fitting the effective potential |T"36j) 
to the pressure and energy density of the pure gauge theory with three colours; 
those are b 3 « 0.9 and bl/[T*db 2 {T*)/dT] « 1, to within 10%. I therefore expect 
that numerically AT C /T* is roughly equal to £q(T~). 

The Eqs. fl2.HI2.12j) seem to indicate that the discontinuity of £q at T c vanishes, 
i.e., that the phase transition turns into a crossover, at a pion mass such that 
6i( m 7r) = 263/27. However, one cannot really extend the 0(bi) estimates to the 
endpoint of the line of first-order transitions because it applies, near T c , only 
if -46 2 (T C ) < &§, which translates into 61 < 63/IO8, see Eq. (I2~TU1) . To find 
the endpoint of the line of first-order transitions I solve numerically for the global 
minimum of (J2.9j) as a function of 62, for given 61, see Fig. 12.21 (left). The numerical 
solution is "exact" and does not assume that 61 is small. I employ 63 = 0.9 to 
properly account for the small latent heat of the pure gauge theory |DP01| lDP02b, 
IDP02al LSD.) OIL IDP02cl IDP02dl ETTTL02] . Also, for b x = 0, this 63 corresponds to 
£ c = 0.6, which is close to the expectation value of the renormalised (fundamental) 
loop for the N c = 3 pure gauge theory [KKPZ021 iPTOil IZan03l IDHL+04j . 
Clearly, for very small 61 the order parameter £q jumps at some b\ = b 2 (T c ), i.e., 
the first-order phase transition persists. (The abscissa is normalised by |&2(T 6 *)| = 
263/9.) I find that the discontinuity vanishes at b\ = 0.026(1), so there is no true 
phase transition for 61 > b\. Nevertheless, I define b c 2 even in the crossover 
regime via the peak of d£o(b 2 )/db 2 . The shift of b c 2 with increasing 61 can now be 
converted into the shift of T c itself by expanding about T*\ 



as already discussed above. I also note that from Fig. 12.21 (left) the susceptibility 
for the Polyakov-loop at its maximum is d£o/db 2 — 3.5, 2, 1.5 for 61 = 0.06, 0.1, 
and 0.126, respectively. That is, the crossover is rather sharp for the values of 61 
shown in the figure. 

Explicit breaking of the Z(3) symmetry of the gauge theory has previously been 
studied in jGK841 IMU951 11+99] . and has been identified as the essential factor 
in determining the endpoint of deconfining phase transitions. Moreover, while 
the term ~ 61 quickly washes out the transition, those studies showed that along 




(2.14) 
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Figure 2.2: Left: The expectation value for the Polyakov-loop, £o(b2(T)), for various 
values of the explicit symmetry breaking coupling, b\. All curves terminate at £o = 
1 44> T = oo. Right: b\ as a function of m n , obtained by matching to three-flavour 
lattice data for T c (m n ). The solid line corresponds to an exponential increase ofb\ with 
decreasing m n , see text. The broken horizontal line displays the endpoint of the line of 
first-order phase transitions in terms ofb\; the intersection with the b\(m w ) curve then 
gives the corresponding pion mass. 



the line of first-order transitions the shift of T c (or, alternatively, of the critical 
coupling /3 C ) is moderate, which agrees with my findings. However, the numerical 
values for the critical "external field" at the endpoint obtained in |M()95[ IA + 99j 
from actual Monte-Carlo simulations can not be compared directly to my estimate 
for b\ because I work here with the renormalised (continuum- limit) loop, not the 
bare loop. Ref. jKLPOlj studied finite-temperature QCD with Nf = 3 flavours 
and various quark masses on the lattice (with improved p4-action), and deter- 
mined the critical (or crossover) temperature as a function of the pion mass. Using 
Eq. IjZHjl one can match A.T C /T* to the data from jKLPOl] to determine hiijn^). 
In other words, I extract the function bx{m % ) required to match the effective La- 
grangian (|2.9jl to T c (m T ) found on the lattice. The result is shown in Fig. !2.2l on 
the right. (Again, the pion mass is normalised to the zero-temperature string ten- 
sion in the pure gauge theory, y/cr ~ 0.425 GeV.) Evidently, the ~ 33% reduction 
of T c from = oo (pure gauge theory) to m^jy/o pa 1 requires only small ex- 
plicit breaking of the Z(3) symmetry for the Polyakov-loop £: I find that b\ < 0.15 
even for the smallest pion masses available on the lattice. This is due to the rather 
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weak first-order phase transition of the pure gauge theory with N c = 3 colours, re- 
flected by the strong dip of the string tension in the confined phase near T~ and of 
the Polyakov-loop screening mass m,£ in the deconfined phase near T c + |KKLL00] ; 
cf . also the discussion in jPOTl IFM ITTPOTI IDP02bl M)P02al M)P02cl IDP02d| . 
Moreover, bi(m n ) appears to follow the expected behaviour ~ exp(—m n ). The ex- 
ponential fit shown by the solid line corresponds to b^m^) = aexp(— b m n /y/a), 
with a = 0.19 and b = 0.47 . Surprisingly, by naive extrapolation one obtains a 
pretty small explicit symmetry breaking even in the chiral limit, b\ ~ 0.2. 
Finally, the endpoint of the line of first-order transitions at b\ = 0.026 (indicated 
by the dashed horizontal line) intersects the curve b^m^) at m n /^a « 4.2. 
For heavier pions the theory exhibits a first-order deconfining phase transition, 
which then turns into a crossover for < 4.2 ^/cr « 1.8 GeV {=m q « 0.9 GeV). 
According to my estimate, the endpoint of the line of first-order transitions occurs 
at AT C /T* w 12%, which is slightly less than a previous (qualitative) estimate of 
26% from ref. jOMOOj . 

Of course, so far my analysis is restricted to pion masses m n /^Jcr > 1. On the 
other hand, one might cross a chiral critical point for some pion mass m^j^Jo < 
1 |OOP94j . Attempting a fit with the model ()2.9j) beyond that point would then 
lead to deviations from b\ ~ exp(— m n ). 



— Ill— 



The improved Hartree-Fock 
approximation 



3.1 Motivation 

In this section I consider the linear cx-model with O(N) symmetry (cf. Sec. 
II. 4j) within the CJT formalism (cf. Sec. 11.5(1 . Of course, it is practically im- 
possible to solve the exact theory (i.e. to include all possible diagrams in the 
effective potential), therefore one has to solve the model within a many-body ap- 
proximation scheme. The most popular among these many-body approximation 
schemes is the Hartree-Fock approximation. In this case, T2 [cf. Eq. (1 1.44)1 ] solely 
contains diagrams of so-called "double-bubble" topology, cf. Fig. 11.41 a-c. Note 
that, in Refs. |Pet991 iLROUl ILRSB001 IRRR03j . this scheme has been termed 
"Hartree" approximation. However, since the exchange contributions with re- 
spect to internal indices are in fact included, it is more appropriate to call it 
"Hartree-Fock" approximation. Neglecting the exchange contributions leads to 
the actual Hartree approximation which, in the case of the O(N) model, has 
also been termed "large- N" approximation |Pet99| ILROOj . In the Hartree-Fock 



approximation, the self-energies, which according to Eq. (jl.50|) are obtained by 
cutting lines in the diagrams for T2, only consist of "tadpole" diagrams, cf. Figs. 
13. li b. c. and !3.2l a.c. For the chiral effective theories, such as the U(Nf) r x U(Nf)t 
and 0(N) models, this approximation scheme has been studied in great detail 
[BTT771 IBK961 IRM981 IXH971 lPet99l iLROOl ILRSB001 IRRR03| . The Hartree-Fock 
approximation is a very simple approximation scheme, since tadpole self-energies 
do not have an imaginary part and, consequently, all particles are stable quasi- 
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particles. Moreover, since tadpoles are independent of energy and momentum, 
the Dyson- Schwinger equations for the propagators reduce to fix-point equations 
for the in-medium masses. 

There are, however, several problems with the Hartree-Fock approximation re- 
lated to the truncation of T 2 . For instance, in the case of the O(N) model, it does 
not correctly reproduce the order of the chiral symmetry restoring phase transi- 
tion and it violates Goldstone's theorem in the sense that the Goldstone bosons do 
not remain massless for nonzero temperatures < T < T c . Several ways to rem- 
edy this shortcoming have been suggested. The simplest one is to neglect sublead- 
ing contributions in 1/N, leading to the Hartree (or large- N, see discussion above) 
approximation |Pet99| iLROOj . Here, the Goldstone bosons remain massless for 
T < T c and the transition is of second-order, as expected from universality class 
arguments |PW84aj . Another possibility to restore Goldstone's theorem is via so- 
called "external" propagators |ABW98I lvHK()2al lvITlTU2cl lvHK02bl IAAB+02| . 
These objects are obtained from the Hartree-Fock propagators by additionally 
resumming all diagrams pertaining to the Random Phase Approximation (RPA), 
with internal lines given by the Hartree-Fock propagators. Another way to re- 
store the second-order nature of the chiral phase transition is to include so-called 
"two-particle point-irreducible" (2PPI) contributions to T 2 jV( ;92ll\7eruT1lBM03| . 
In this work, I am not concerned with these formal shortcomings of the Hartree- 
Fock approximation: I focus exclusively on the case realized in nature where chiral 
symmetry is already explicitly broken by (small) nonzero quark masses, such that 
the pion is no longer a Goldstone boson and assumes a mass m n ~ 139.5 MeV 
[EH O + 04j . My goal in this work is to include the nonzero decay width of the 
particles in a self-consistent many-body approximation scheme. To this end, one 
has to go beyond the Hartree-Fock approximation and add other diagrams to T 2 
which, upon cutting lines according to Eq. (|1.5()jl . yield self-energies with nonzero 
imaginary part. The most simple of such diagrams, and the ones considered in 
the following, are those of the so-called "sunset" topology, cf. Fig. ll.4l d.e. leading 
to the self-energy diagrams shown in Figs. 13.11 d,e and 13.21 d. 
The self-energies arising from the sunset diagrams in T 2 depend on the energy and 
the momentum of the incoming particle. Therefore, the Dyson- Schwinger equa- 
tions no longer reduce to fix-point equations for the in-medium masses. Since the 
self-energies now have a nonzero imaginary part, implying a finite lifetime of the 
corresponding particles, the spectral densities are no longer delta-functions. It is 
therefore convenient to rewrite the Dyson-Schwinger equations for the propaga- 
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tors into a set of self-consistent integral equations for the spectral densities which 
has to be solved numerically. 



3.2 The Dyson-Schwinger and the condensate 
equations 

In this section I discuss the application of the improved Hartree-Fock approxi- 
mation to the O(N) model. The linear a model with O(N) symmetry and its 
effective potential is already introduced in Chapter 11.41 As mentioned above, in 
the standard Hartree-Fock approximation one takes into account only the three 
possible double-bubble diagrams, V 2 , V 2 , and V 2 of Eqs. 11.291 in the effective po- 
tential. As explained, in order to include the nonzero decay width of the particles 
one has to go beyond the standard Hartree-Fock approximation by additionally 
including the two sunset diagrams, V 2 and V 2 of Eqs. 11.301 The complete ex- 
pression for V 2 is obtained by the sum of all contributions, 



v 2 = v 2 a + v 2 b + V 2 C + v 2 d + v 2 e 



(3.1) 



The expectation values of the one- and two-point functions in the absence of 
external sources, a and S, V, are determined from the stationary points of V, 
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5=a,S=S,P=V 
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a=a,S=S,P=V 

leading to an equation for the scalar condensate a, 

4A 4A r 
H = ^a+-a 3 + -a J [35(g) + (N - 1) V{Q)\ 



, (3.2) 



a=a,S=S,P=V 
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(3.3a) 



a=a,S=S,P=V 



where 



5V, 
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a=a,S=S,P=V 



'4V 
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a 



2(N-1) / / S(L)V(Q)V(L + Q) 



+ 3! / / S(L)S(Q)S(L + Q) 

JL JQ 



(3.3b) 



and to the Dyson-Schwinger equations for the scalar and pseudoscalar propaga- 
tors, 



S-\K-a) = S- 1 (K;a) + Z[K;a] , 
V-\K-a) = p-\K;a) + U[K;a] . 



(3.3c) 
(3.3d) 
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Figure 3.1: The self-energy of the a-meson. The diagrams b and c are tadpole contri- 
butions generated by cutting an internal cr-line in the double-bubble diagrams b and c 
in Fig. 11.41 The diagrams d and e are one-loop contributions generated by cutting an 
internal a-meson line in the sunset diagrams d and e of Fig. \1.4l Lines and vertices as 
in Fig. EH 



Here I introduced the self-energies of the scalar and pseudoscalar fields, 



E = S a + S fe + S c + S d + S e , II = IT + IT + IT + U d + IT 



(3.4) 
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a c d 




Figure 3.2: The self-energy of the pion. The diagrams a and c are tadpole contributions 
generated by cutting an internal pion line in the double-bubble diagrams a and c in 
Fig. 11.41 The diagram d is the one-loop contribution generated by cutting an internal 
pion line in the sunset diagram d in Fig. 11.41 Lines and vertices as in Fig. 11.41 
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The calculation of the self-energy contributions Eqs. (|3.5|) and ()3.6|) corresponds 
to opening internal lines in the diagrams of Fig. 11.41 Employing this procedure 
to the double-bubble diagrams leads to the tadpole contributions for the self- 
energies of the cx-meson (see Figs. 13.11 b. c) and the pion (see Figs. 13.21 a. c). 
This defines the standard Hartree-Fock approximation. Additionally, the sunset 
diagrams of Figs. I1.4| d. e lead to the one- loop contributions shown in Figs. I3.l| d. 
e and !3.2l d. The latter contributions depend on the energy and the momentum of 
the particles and lead to nonvanishing imaginary parts for the self-energies. The 
explicit calculation of the self-energies (|3.5j) and (|3.fij) is discussed in Appendix 
A.l. The next step is to rewrite the set of Dyson-Schwinger equations (|3.3c|) and 
(j3.3d|) in terms of a set of integral equations for the spectral densities. In general, 
the spectral densities for the a-meson and the pion are defined as 

p a (u,k) = 2 Im S(u + irj, k; a) , (3.7a) 
p n (co,k) = 2 Im P(lu + irj, k; a) . (3.7b) 

If the imaginary parts of the self-energies are zero, as in the standard Hartree-Fock 
approximation, all particles are stable quasiparticles, i.e., their spectral densities 
become delta- functions with support on the quasiparticle mass shell, 

p ff (w,k) = 27rZ CT [u; CT (k),k] {8[u - u a {k)) - 5[u + u a {k)}} , (3.8a) 
PlT (uj,k) = 2 7 rZ 7r K(k),k] {5[lu - luM] - 5[lu + luM}} , (3.8b) 

where a> CT)7r (k) are the quasiparticle energies for a-meson and pions, defined by 
the positive solutions of 

uj 2 a (k) = k 2 + ml(a) + ReJ:[uj a (k),k;cr} , (3.9a) 
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and 
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(3.9b) 

(3.10a) 
(3.10b) 



are the wave-function renormalization factors on the quasi-particle mass shell. 
Since the real parts of the self-energies are even functions of energy, these factors 
are also even in c<j CTj7r (k). In the standard Hartree-Fock approximation, the (real 
parts of the) self-energies do not depend on K^, such that 



where 



M 2 (a) = ntL(a) + ReE(cr) , M 2 (a)=m 2 (a) + ReU(a) 



are the effective cr-meson and pion masses, and 

1 



Z a [u a {k)} 



2u a (k) 



Z^u^k)} 



2u w (k) 



(3.11) 



(3.12) 



(3.13) 



For nonvanishing imaginary parts of the self-energies the spectral densities assume 
the following form: 



2 Im £(co>, k; a) 



[u 2 -k 2 — m 2 (a) - Re k; a)} 2 + [Im E(u, k; a)] s 
2Imn(w,k;a) 



, (3.14a) 



- - k 2 - m 2 n (a) - Re U(cu, k; a)] 2 + [lmU(uj, k; a)] 2 - (3 ' 14b) 



The general calculation of the imaginary parts of the self-energies is discussed 
in Appendix A.l. It turns out that they do not depend on the direction of 3- 
momentum k, but only on the modulus k. The tadpole diagrams of Figs. 13. li b, 
c (for the cx-meson) and Figs. 13.21 a, c (for the pion) do not have an imaginary 
part. Thus [see Appendix, Eq. (|A.18|) ]. 
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x 



roo 

/ dqi qi dq 2 q 2 © (|<?i - 92 1 < k < q x + q 2 ) 
Jo 

[2(N - 1) pn(ui, q 1 )p 1T {u) 2: 92) + 3 ■ 3! p a (ui, gi)p CT (cj 2 , 92)] ,(3. 15a) 
Imn(w,k;a) = Imn d (w,k;(r) 

1 1 Z" 00 

2(27rj3 fc / duJl dU2 ^ + + ^ U2 ^ u ~ Ul ~ U2 > 

roo 

/ dq 1 q 1 dq 2 q 2 Q (\q x - q 2 \ < k < q 1 + q 2 ) 
Jo 

X Ap (J (u 1 ,q 1 )p 1T (u 2 ,q 2 ) , (3.15b) 

where f(uS) = [exp(u/T) — 1] _1 is the Bose-Einstein distribution function. These 
expressions are not ultraviolet divergent and thus do not need to be renormalised. 
For the real parts, I employ the following approximation. One only takes into ac- 
count the Hartree-Fock contributions, arising from the tadpole diagrams Figs. 13.11 
b, c (for the cr-meson) and Figs. 13.21 a. c (for the pion). As discussed above, these 
contributions do not depend on energy and momentum. The real parts of the 
diagrams in Figs. 13.11 d, e and 13.21 d are functions of energy and momentum and 
are much harder to compute, involving an additional numerical integration as 
compared to the respective imaginary parts ()3.15j) . Therefore, in the present first 
step where I focus exclusively on the effects of a nonzero decay width, I neglect 
them. Consequently [see Appendix, Eq. (jA.9|) ]. 

ReE(a) = Re S b (a) + Re S c (a) 



4A 4 



N (2tt) 3 

Rell(cx) = ReIT(a) + ReIT( ( T 
4A 4 



dudqq f(u) [3 p a (u, q) + (N - 1) p^(u, q)] , (3.16a) 



duodqq 2 f(uj) [p <7 (u,q) + (N + l)p^u,q)\ . (3.16b) 



N (2tt) 3 

Since the real parts do not depend on energy and momentum, they only modify 
the masses of the a-meson and the pion, cf. Eq. ()3.12j) . Note that I only consider 
the temperature-dependent contributions in Eqs. ()3.16|) . The omitted vacuum 
parts are ultraviolet divergent and have to be renormalised. The proper way to 
perform renormalization within the CJT-formalism was first presented in Refs. 
[vHK02a, vHK02c t lvHK02b"] . In the standard Hartree-Fock approximation for 
the O(N) model, this leads to the same equations that one obtains applying the 
renormalization procedure of Ref. [LROOJ. In that paper, the dependence of the 
results on the renormalization scale was studied in detail. It was found that sim- 
ply neglecting the vacuum parts instead of applying the proper renormalization 
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procedure leads only to quantitative, but not to qualitative changes in the results. 
Since I am interested in the effect of nonzero imaginary parts for the self-energies 
and not in issues of renormalization, I simply neglect the vacuum parts. This 
prescription is also used to compute the tadpole and the sunset contributions in 
Eq. (J3.3a|) for the condensate. 

After obtaining real and imaginary parts (|3.16|) and (|3.15|h respectively, one in- 
serts them in the expressions (|3.8|) (if the imaginary part is zero) or ()3.14|) (if 
the imaginary part is nonzero) for the spectral densities. These spectral den- 
sities can then be used to again evaluate the real and imaginary parts of the 
self-energies. This defines an iterative scheme to self-consistently solve for the 
spectral densities as a function of energy and momentum. A convenient starting 
point for this scheme is the standard Hartree-Fock approximation, i.e., neglecting 
the imaginary parts altogether. 

The self-consistently computed spectral densities obey a sum rule jTBOO , 



In my calculation, there are two reasons why this sum rule may be violated. One 
is because I neglected the real parts of the self-energy diagrams in Figs. 13.11 d, e 
and 13.21 d. The other is due to the numerical realization of the above described 
iterative scheme. Numerically, one has to solve for the spectral density on a 
finite, discretized grid in energy-momentum space. If the imaginary part of the 
self-energy becomes very small, the spectral density converges towards a delta 
function. The support of the delta function may be located between grid sites, 
which causes a loss of spectral strength, which in turn violates the sum rule. 
My prescription to restore the validity of the sum rule is the following. One 
first checks whether the imaginary part is smaller than the grid spacing (in my 
calculations, 10 MeV in both energy and momentum direction) at the location 
of the quasiparticle mass shell, ou = ou a>7r (k). If this is the case, one adds a 
sufficiently wide numerical realization of the delta function, <5 num , to the original 
spectral density, p(u, k) — > p(ui, k) + c ■ 5 num [u — LU aj7r (k)], where c is a constant 
that is adjusted so that the sum rule is fulfilled on our energy- momentum grid, 



where oo max is the maximum energy on the energy- momentum grid. 

On the other hand, if the imaginary part turns out to be sufficiently large, I 

presume that a possible violation of the sum rule 1)3.18)1 arises from neglecting 




(3.17) 




(3.18) 
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the real parts of the self-energy diagrams in Figs. 13. li d, e and I3.2l d. In this case, 
one multiplies the spectral density by a constant, p — > d ■ p, where d is adjusted 
so that the sum rule ()3.18|) is fulfilled. (By comparing the results to the case 
d — 1, I found that this somewhat ad hoc correction procedure does not lead to 
major quantitative changes.) 

Restricting the sum rule to a finite range in energy as in Eq. (J3.18|) . however, 
causes the following problem. If the decay width of the particles is very large and 
consequently the spectral density a rather broad distribution around the quasi- 
particle mass shell, there will be parts which lie outside the energy-momentum 
grid. I could estimate the magnitude of this physical effect if I knew the behaviour 
of the spectral density at energies uj > u max . This is possible at zero temperature 
(with the help of Weinberg's sum rules), but not at nonzero temperature. Here, I 
simply assume that this effect is sufficiently small to be neglected, i.e., I assume 
that a possible violation of the sum rule (|3.18|) is due to the two above mentioned 
artifacts and accordingly perform the correction of the spectral density. 
Another physical effect which causes a loss of spectral strength is if the quasipar- 
ticle energy ui a ^{k) > u milx . In this case, I do not perform the correction of the 
spectral density as described above. This occurs at large momenta k, close to the 
edge of the energy-momentum grid. I checked that, in the numerical calculation 
of the integrals in Eqs. (j3.15|) . the integrands become sufficiently small in this 
region, so that the imaginary parts are not sensitive to this effect. 
This concludes the discussion of the improvement of the standard Hartree-Fock 
approximation by self-consistently including nonzero decay widths. 

3.3 Results 

In this section I present the numerical results for the 0(4) linear sigma model 
obtained in the Hartree-Fock approximation improved by including nonzero de- 
cay widths, as discussed in Sec. 3.2. The condensate a is shown in the left 
part of Fig. 13.31 as a function of temperature for the standard and the improved 
Hartree-Fock approximation. The qualitative behaviour is similar in the two 
approximations: the condensate drops significantly with temperature indicating 
the restoration of chiral symmetry. Since I consider the case of explicit 0(4) 
symmetry breaking by taking H ^ in Eq. (jl.22)l . the chiral phase transition 
is a crossover transition. Nevertheless, one can define a transition temperature, 
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Figure 3.3: The values for the condensate (left panel) and the effective masses of 
the a-meson and the pion, (right panel) as functions of temperature. The values are 
calculated in the standard Hartree-Fock approximation (dashed and solid lines) and in 
the improved Hartree-Fock approximation (symbols) as discussed in Sec. 3.2. 



T x , as the temperature where the chiral susceptibility da/dT assumes a maxi- 
mum. Quantitatively, the inclusion of a nonzero decay width lowers T x by about 
20% as compared to the standard Hartree-Fock approximation. The transition 
temperature T x ~ 175 MeV agrees within error bars with recent lattice results, 
T x ~ 172 ±5 MeV for the two- flavour case |LP03j . (Note, however, that the latter 
value is extracted from an extrapolation to the chiral limit, while my results are 
for the case of explicit symmetry breaking, i.e., for nonzero quark masses.) 

For the solution of the condensate equation (|3.3ajl the relative magnitude of the 
contribution from the sunset diagrams, Eq. (j3.3b|) . is negligibly small, of order 
~ 10~ 4 , and thus can be safely neglected. In turn, not having to compute the 
integrals pertaining to the double loop, cf. Eq. (jA.24|) . considerably speeds up 
the computation. 

In the right part of Fig. 13.31 I show the effective masses M a = m 2 a + Re £ and 
M n = m 2 + Rell of the cr-meson and the pion as functions of temperature in 
both approximation schemes. Since the decay width of the pion remains com- 
paratively small, cf. Fig. 13.41 the mass of the pion does not change appreciably 
when taking the nonzero decay width into account. On the other hand, the large 
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Figure 3.4: The decay width of the a-meson, T a = ImE(u a )/u a , and the pion, 
T-X = ImU(uj n ) / uj n , in the improved Hartree-Fock approximation as a function of the 
temperature at momentum k = 325 MeV. 



decay width of the cr- meson at temperatures below T x , cf. Fig. 13.41 does influence 
the mass: in this range of temperatures the mass exhibits a stronger decrease 
with temperature in the improved Hartree-Fock approximation. At large tem- 
peratures T > T x , the decay width of the a-meson becomes negligibly small, cf. 
Fig. 13 A\ and the mass approaches the value computed in the standard Hartree- 
Fock approximation. Both a-meson and pion masses approach each other above 
T x , indicating the restoration of chiral symmetry. 

In Fig. 13.41 1 show the decay widths of a-mesons and pions, defined as |Wel83} 

iLBnnl 

ImEMM 1 M ) r^E^flM, (3.19) 

where uj a ^{k) is the energy on the quasi-particle mass shell, cf. Eq. 1)3.9}) . At 
small temperatures, due to the possible decay of a a-meson into two pions, the 
decay width of the former is large, of the order of its mass. Note that the 
value T a ~ 460 MeV obtained here at T = is completely determined by the 
parameters m 2 ,A, and H of the 0(4) model, i.e., without adjusting any addi- 
tional parameter. It is reasonably close to the experimentally measured value 
r CT ~ (600 - 1000) MeV )EHO+04j . The a decay width increases with tempera- 
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ture up to a maximum of about 500 MeV at a temperature T ~ 100 MeV and 
then decreases rapidly. The decay width of the pion vanishes at T = 0. It in- 
creases at nonzero temperature and assumes a maximum at about T ~ 180 MeV. 
At this temperature the decay width is about half of the corresponding pion mass. 
It decreases rapidly at higher temperature. Although the decay widths of both 
particles decrease at high temperatures, they do not become degenerate, the de- 
cay width of the cr-meson remains about a factor of 8 larger than that of the pion. 
This difference can be traced to the symmetry factors multiplying the one-loop 
self-energies for the a-meson and pion, cf. Eqs. (j3.15|) . For the cr-meson, there 
is a factor 2(N — 1) in front of the pion loop and a factor 3 • 3! in front of the 
cx-loop, so that the overall factor is ~ 2(N — 1) + 3 ■ 3! = 24. For the pion, there is 
only a mixed cr-pion loop with a symmetry factor 4. From this simple argument 
one already expects that the decay width of the cr-meson is about a factor of 6 
larger than that of the pion. The remaining difference comes from the fact that 
the self-consistently computed spectral densities of cr-meson and pion under the 
integrals in Eqs. (|3.15|) are also different, cf. Fig. 13.71 

One might argue that, at asymptotically large temperatures, effects from chiral 
symmetry breaking can no longer play a role and the decay widths, as well as the 
spectral densities, of cx-meson and pion should become degenerate. This is true 
in the chiral limit, where there is a thermodynamic phase transition between the 
phases of broken and restored chiral symmetry, and where a = in the latter 
phase. Here, however, I consider the case of explicitly broken chiral symmetry, 
where a > 0, even when the temperature is very large. Since the decay widths are 
proportional to a 2 , they also do not vanish at large temperature. The difference 
in the symmetry factors for the self-energies of cr-mesons and pions then leads to 
different values for the decay widths and spectral densities. 
The self-consistently calculated spectral densities of the cr-meson and the pion 
as functions of the external energy uj and momentum k are shown in Fig. 13.51 
and Fig. 13.61 for different temperatures, T = 80, 160, 240, and 320 MeV. For a 
detailed discussion, let us fix the momentum at k = 325 MeV and consider the 
spectral densities as functions of energy u for different temperatures, cf. Fig. 13.71 
At all temperatures, the pion spectral density (dotted line) exhibits a pronounced 
peak on the mass shell, at u n (k) = \Jk 2 + M 2 (a). When the temperature is below 
200 MeV, such that M w ~ m n = 139.5 MeV, cf. Fig. EE the peak is located 
at 0^(325 MeV) ~ 350 MeV. Above T ~ 200 MeV, M n increases significantly 
with temperature, and the position of the peak shifts towards larger energies, 




Figure 3.5: The spectral density of the cr-meson as a function of energy uj and momen- 
tum k at temperatures 80, 160, 240, and 320 MeV. 




Figure 3.6: The spectral density of the pion as a function of energy uj and momentum 
k at temperatures 80, 160, 240, and 320 MeV. 
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Figure 3. 7: The spectral density of the cr-meson and pion as a function of energy lo at 
temperatures 80, 160, 240, and 320 MeV. The momentum is k = 325 MeV. 
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a; 7r (325 MeV) ~ 500 MeV. The broadening of the peak is due to scattering of the 
pion off (T-mesons in the medium. 

In contrast to the pion spectral density, for temperatures below ~ 170 MeV the 
a spectral density (full line) does not exhibit a peak at the mass shell, u a (k) = 



yk 2 + M 2 (a). The reason is that u a (k) is still sufficiently large to allow for the 
decay into two pions. Consequently, in this temperature range the a spectral 
density is very broad. On the other hand, for temperatures above ~ 170 MeV, 
where ou a (k) drops below 2M n , the two-pion decay channel is closed and the a 
spectral density develops a distinct peak, whose width is due to scattering of the 
cr-mesons off pions and other a-mesons in the medium. 

Two other features of the spectral densities shown in Fig. 13.71 are noteworthy. The 
first is the region below the light-cone, K 2 = uj 2 — k 2 < 0, where the mesons are 
Landau-damped. The second is the two-particle decay threshold. For a — > rnr, 
this threshold is located at uj ~ 2M n , for a — > era, it is at u> ~ 2M„, and 
for 7i — > cm at uj ~ M a + M n . The threshold is most clearly seen at large 
temperatures, e.g. T = 320 MeV, when both particles become degenerate in 
mass, M a ~ M w ~ 400 MeV, and the threshold is at u ~ 900 MeV. 




-IV- 



The improved Hartree 
approximation 

4.1 Motivation 

In the standard Hartree (or Hartree- Fock) approximation of the O(N) model, 
one only takes into account the double-bubble diagrams. In these approximations 
the self-energies of the particles are only real valued, therefore no decay width 
effects are included. The difference between the Hartree and the Hartree-Fock 
approximation is, that in the Hartree approximation all terms of order ~ 1/N 
are neglected on the level of the Dyson- Schwinger and the condensate equations. 
In the last chapter I presented the so-called improved Hartree-Fock approxima- 
tion, which takes additionally into account sunset type diagrams. This leads to 
4-momentum dependent real and imaginary parts in the Dyson-Schwinger equa- 
tions. Indeed, in Chapter III, I neglected the 4-momentum dependent real parts 
of the Dyson-Schwinger equations for simplicity. In the present Chapter, I include 
them, although in the Hartree approximation |Rod05j . This leads to a vanishing 
imaginary part of the pion self-energy (because it is of order ~ 1/A^"), i.e., the 
spectral density of the pion is just a delta function. As shown in chapter III, the 
inclusion of the pion decay width does not lead to a major change of the results, 
therefore this is still a good approximation for the pion. The contributions from 
the a — > 2a decay, in the self-energy of the cr-meson, is also of order ~ 1/N, and 
vanishes, but (a part of) the contribution from the a — > 2tc decay remains. There- 
fore, its spectral density has a nonzero width, as expected for the cx-meson with 
a very large (vacuum) decay width, T a = (600 — 1000) MeV. In the following, I 
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(J -lilcbUIl 

vacuum mass 


explicit chiral symmetry breaking 


chiral limit 


= 139.5 MeV, f v = 92 A MeV 


= MeV, / ff = 90 MeV 


H = (121.6 MeV) 3 


H = 


m a = 400 MeV 


A = 8.230 


A = 9.877 




u 2 = -(225.41 MeV) 2 


/x 2 = -(282.84 MeV) 2 


m a = 600 MeV 


A = 19.043 


A = 22.222 




^ = -(388.34 MeV) 2 


u 2 = -(424.264 MeV) 2 


m a = 800 MeV 


A = 36.341 


A = 39.506 




^ = -(539.27 MeV) 2 


a 2 = -(565.685 MeV) 2 



Table 4.1: The masses and decay constants at vanishing temperature and the corre- 
sponding parameter sets for the 0(4) linear sigma model for the two symmetry breaking 
patterns studied here. 

call this the improved Hartree approximation. 

As discussed in Sec. 1.4, the parameter H is given as a function of the vacuum 
mass m n , and the vacuum decay constant f n of the pion : H = m%f n . In this 
chapter I compare the chiral limit, = and f n = 90 MeV, with the case of 
explicit chiral symmetry breaking, m n = 139.5 MeV and f n = 92.4 MeV. The 
bare mass n 2 and the four-point coupling A depends additionally on the vacuum 
mass of the a-meson m G : fi 2 = —(m 2 — 3m 2 )/2, and A = N(m 2 — m 2 )/(8f 2 ). 
As mentioned above, the decay width of the er-meson in vacuum is very large, 
therefore its vacuum mass is not well defined, m a = (400-1200) MeV EHO+04. 
I compare the results for m a = 400, 600, and 800 MeV. The parameter sets of 
the model, for N = 4, are summarised in table 14.11 

4.2 The Dyson-Schwinger and the condensate 
equations 

In the following I repeat the derivation of the condensate and the Dyson- 
Schwinger equations in the CJT formalism, as discussed in chapter III, and show 
how the equations simplify in the improved Hartree approximation. The dia- 
grams included in the effective potential are the three double-bubble diagrams 
[~ (/ P) 2 , ~ (/ S) 2 , and ~ / P j S] shown in Figs. I1.4| a. b, c, and the two sunset 
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diagrams (~ / / SPP and ~ / / SSS) shown in Figs. 11.41 d, e. To get the ex- 
pectation values for the one- and two-point functions in the absence of external 
sources, a, S, and V, one has to find the stationary points of the effective poten- 
tial ()1.25)) . Minimisation of the effective potential with respect to the one-point 
function 

AT/ 

0, (4.1) 



5V_ 
IF 



a=a,S=S,P=V 

leads to an equation for the scalar condensate a, 

4A , 4A 



H = fi 2 a + -a 3 + ^aJ Q [3S(Q) + (N-l)V(Q)] 

2(N-1) f f S(L)V(Q)V(L + Q) 
Jl jq 



4A 



3! / / S(L)S(Q)S(L + Q) 

JL JO 



(4.2) 



Using the fact that a 2 ~ iV [cf. Eq. (jl.26j) ]. this equation becomes in the large- N 
limit, 



a 2 + N [ V{Q) 

JQ 



(4.3) 



The minimisation of the effective potential with respect to the two-point func- 
tions, 



sv_ 
Is 



5V 



= 0, jp =0, (4.4) 

a=a,S=S,P=V a=a,S=S,P=V 

leads to the Dyson-Schwinger equations for the scalar and pseudoscalar propaga- 
tors, S and Q, 



S-\K;a) = S-\K-a) + Y,{K-a) , 

Here I introduced the self-energy of the scalar 

SV 2 



(4.5a) 
(4.5b) 



E(K;a) = 2 



5S 



a=a,S=S,P=P 



4A 



3 / S(Q) +(N-1) [ V(Q) 

JQ JQ 



+ 



N 



2(N 



V{K-Q)V{Q) 



+3-3! [ S(K-Q)S(Q) 

JQ 



(4.6a) 
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and pseudoscalar fields 

2 5V 2 



N- 1 5P 
4A 

iV 



a=cr,S=S,P=V 



j Q S(Q)UN + l)J Q V { Q) 



In the large- limit, all a-meson tadpoles and a part of the pion tadpole (~ / S 
and ~ / V) vanish, and only a part of the pion term (~ / VV) remains, 

V(K;a) = j Q V{Q) +(^j2N j Q V{K-Q)V{Q), (4.7a) 

U(K;a) = J q V{Q). (4.7b) 
The tadpole contributions have no imaginary parts, therefore 

(A \ \ 2 

2iVIm J^V(K-Q)V(Q) , (4.8a) 

Imn = 0. (4.8b) 



The imaginary part of the a-meson self-energy depends on the 4-momentum 
vector, K = (k, u), but the real part can be split into terms which do not depend 
on K, [Re£]i, arising from the tadpole contribution, and terms which are 4- 
momentum dependent, [Re S(K, <r)] 2 , arising from the sunset diagram, 

ReS(K ;( j) = [Re£]i + [Re£(K;<7)] 2 , (4.9) 

where 

4A r 

[ReE]! = —N J V(Q), (4.10a) 

[ReE(K;a)] 2 = (^pj 2 NReJ^ V(K - Q)V(Q). (4.10b) 

As mentioned above, the 4-momentum dependent terms in the pion self-energy 
vanish in the large- N limit. The 4-momentum independent term of the real part 
of the pion self-energy is the same as for the cr-meson, 

4A r 

Ren=[ReE] 1 = — N J V{Q) . (4.11) 
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I want to calculate the spectral density of the cr-meson, p a . To this aim, I rewrite 
the Dyson-Schwinger equations (|4.5J) and the equation for the chiral condensate 
(14. 3(1 as functions of the spectral density of the pion p n . In the large- N limit, the 
imaginary part of the pion vanishes, and therefore the spectral density is just a 
delta function 

p v (u, k) = — ^- {5[u - u^k)) - 5[u + u v (k)]} , (4.12) 
with support on the quasiparticle energy for the pion u n (k), 



u,(k) = ^k 2 + M 2 (a), (4.13) 

where I defined an effective mass for the pion M v , 

M^(cr) = ml(a)+ReU. (4.14) 

In the chirally broken phase (a ^ 0) the imaginary part of the cr-meson is nonzero, 
therefore the spectral density assumes the following form: 

2 Im £(u;, k; a) 

Pff ^' 3 = ~ [u 2 - k 2 - mKa) - Re Z(u, k; cr)] 2 + [ImS(w, k; cr)] 2 ' (4 ' 15) 

In the chirally restored phase (cr = 0) the 4-momentum dependent parts of the cr- 
meson self-energy vanish, hence also the spectral density of the cr-meson becomes 
a delta function, 

p a (u, k) = {5[uj - u a {k)} - 5[u + u a (k)]} , (4.16) 
uj a [k) 

with support on the quasiparticle energy for the cr-meson u a (k), 

tu a (k) = ^fk 2 + M 2 (a), (4.17) 

where I defined an effective mass for the cr-meson M a , 

M*(<7)=7r£(<7) + [Re£]i. (4.18) 

Note that in the chirally broken phase (cr ^ 0), the quasiparticle energy for the 
cr-meson uj a (k) is given by the solution of 

ul(k) - k 2 - M 2 a {a) - [ReSK(k),k; a)] 2 = 0. (4.19) 

The imaginary and the real parts of the self-energies, Eqs. ([4.8a]) . (j4.10b|) . and 
f)4.11j) . can be written as functions of the pion spectral density. Note that for 
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the real parts, I only consider temperature-dependent contributions and neglect 
the (ultraviolet divergent) vacuum parts, which is a simple way to renormalize 
the integrals. The imaginary parts are not ultraviolet divergent and thus do not 
need to be renormalised. The 4-momentum independent term is just the standard 
tadpole integral, cf. Eq. (|A.26j) . 

ReIT= [ReE]! = — / dq ^M?)] -1 ^?)]. (4.20) 

Z7T Z JO 

The derivation of the equations for the 4-momentum dependent terms is given in 
appendix A. 2, 

/4Acx\ 2 N 1 r°° 

ImS(cu,k) = (— J — -j Q dqq[u- K {q)]- l Q{\q Q -q\ < k < q + q) 

x {1 + /K(g )] + /K(g)]}, (4.21a) 
[Re£(u;,k)] 2 = ^— j — -P J q dq iqi dq 2 q 2 

x © (|?i - fel < k < qi + q 2 ) [^(gi)^^)] -1 

x [ /K(gi)] + /K(g 2 )] , + (421b) 

, + /K(g 2 )] + -/KOzO] + /K(g 2 )] 1 



^(gi) - LO^{q 2 ) - uj Wn(qi) ~ ^{q 2 ) + u J 
where /(u;) = l/[exp(u/T) — 1] is the Bose-Einstein distribution function, q 



[uj — ^^(g)] 2 — M 2 , and P / . . . denotes the principal value of the integral. Note 
that 6 (|gi — q 2 \ < k < gi + q 2 ) /k = 25 (qi — q 2 ) in the limit k — > 0, which can be 
used to perform the g-integration, cf. Eqs. ()A.38|) . and (jA.44|) . An appropriate 
way to perform the principal value numerically is discussed in the appendix, cf. 
Eqs. flA~l3j) . and 

The spectral densities have to obey a sum rule [LB00J, 

00 da; 

— u Pa>n {uj,k) = l. (4.22) 

-oo Z7T 

In the improved Hartree approximation, the spectral density of the pion (|4.12|) 
is just a delta function, and normalised in such way that this sum rule is a 
priori fulfilled. As mentioned above, this is not the case for the cx-meson spectral 
density in the chirally broken phase. The main reason for a possible violation of 
the sum rule is due to neglecting terms of the order ~ 1/N in the self-energy, 
which leads to a loss of spectral strength. I found that the inclusion of the 
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4-momentum dependent real part [Re£]2 is close to negligible for the validity 
of the sum rule. Other possible problems arise from the numerical realisation 
of the spectral density on a finite uj max x k max energy-momentum grid, where 
Umax > M and k max > \k\ are the boundaries of the grid. If the imaginary part 
of the self-energy becomes very small, the spectral density becomes more and 
more a delta function, which one has to realise numerically. On the other side, if 
this imaginary part is very large, i.e., the width of the spectral density becomes 
very large, one has to use a large u> max , otherwise one loses too much spectral 
strength for energies uj > u) max . To minimise these numerical problems, I use a 
rather large and fine (quadratic) energy-momentum grid, with u) max = k max = 2 
GeV and a lattice spacing of 5 MeV. 

As discussed in Chapter III, if the sum rule is not fulfilled, I use the following 
way to restore it. If the imaginary part of the cr-meson self-energy is very small, 
I add a numerical realisation of a delta function 5 nmn to the spectral density: 
p' a {uJ, k) — > Paiyj, k) +ci • 5 nU m[^ — w a (kj\. If the imaginary part of the self-energy 
is large enough (compared with the lattice spacing) I just multiply the spectral 
density by a factor: p' a (u, k) — > 02- p a {oo, k). The constants c% and C2 are adjusted 
in a way that p' a fulfills the sumrule on the energy-momentum grid 



The numerically scheme for the improved Hartree approximation is the following. 
At first, one solves the standard Hartree approximation, i.e., the condensate and 
the Dyson-Schwinger equations, Eqs. (|4.3|) and (|4.5|) with [ReS] 2 = [ImS] 2 = 0, 
to obtain the chiral condensate and the effective mass of the pion, a and M n . With 
these results, one calculates the 4-momentum dependent real and imaginary parts 
(j4.21|) of the cr-meson self-energy, and finally the spectral density, Eq. ()4.15|) or 
(j4.16j) . The decay width of the cr-meson T CT is defined as jWe!83l [LB00J 




(4.23) 



F*(k) 



ImS[cL> .(k),k; a\ 
uj a (k) 



(4.24) 



4.3 Results 



In this section I present the numerical results for the linear sigma model with 
0(4) symmetry in the improved Hartree approximation as discussed in the last 
section, for the parameter sets given in table I4~T1 
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Figure 4.1: The chiral condensate a (left column) and the effective masses of the o- 
meson and the pion, M a and M n (right column), as functions of temperature T and 
a-meson vacuum mass m a . In the upper row the results with explicit chiral symmetry 
breaking are shown and in the lower row the results in the chiral limit. 



The mass 

I start the discussion of the results with the chiral condensate a, given by the 
solution of Eq. (j4.3|) . and the 4- momentum independent effective masses of the 
ex-meson and the pion, M a and M n , as defined in Eqs. (j4.18J) and (j4.14|) . Note, 
however, that (in the chirally broken phase, <j^0) the mass of the cx-meson (the 
energy where the a-meson spectral density assumes its maximum) is additionally 
modified by the 4-momentum dependent real part of the cx-meson self-energy 
[Re£(K; cx)]2 in Eq. (j4.10b|) . which is discussed later. 

In the upper row of Fig. 14.11 the results for the chiral condensate a (left column) 
and the effective mass for the cx-meson and the pion, M a and (right column), 
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Figure 4.2: The 4-momentum dependent real part of the a-meson self-energy [Re £]2, in 
the energy-momentum plane, for the case with explicit chiral symmetry breaking (left ) 
and the corresponding result for the chiral limit (right) at a temperature of T = 120 
MeV. 



are shown as functions of T and m CT , in the case with explicit chiral symmetry 
breaking. The results show the behaviour of a crossover transition. Neither the 
chiral condensate nor the mass of the scalar particle M a become equal to zero, 
even for high temperatures. Thus, M a and M n become only approximatively 
degenerate for large temperatures. In the lower row of Fig. I4.1l the corresponding 
results for the chiral limit are presented. The results show the behaviour of a 
second-order phase transition. The condensate a and the mass of the scalar 
particle M a becomes equal to zero at a critical temperature T x , therefore M a = 
M w for T > T x . The condensate (accordingly T x ) does not depend on the vacuum 
mass 77V. This can be understood as a consequence of the condensate equation 
flOJ) in the chiral limit (H = 0), 



2 4A 
" + N 



a 2 + N [ V(Q) 

JQ 



(4.25) 



For m n = the integral can be performed analytically LR00J, 



m 2 4A T 2 
= _I^ + 1V + 4A— . (4.26) 
2 N 12 K J 
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Figure 4.3: The 4-momentum dependent real part of the o-meson self-energy [ReX] 2 , 
at a hxed momentum ofk = l GeV, for the case with explicit chiral symmetry breaking 
(full line) and the corresponding result for the chiral limit (dotted line) at a temperature 
ofT = 120 MeV. 

Using (jl.28aj) . this leads to an (m^-independent) equation for the condensate: 



that determines the critical temperature (where a = 0) to be T x = yl2/Nf v . 
This agrees with the numerical results for the 0(4) model, T x = v^/tt ~ 160 MeV. 
The 4-momentum dependent real part of the ex- meson self-energy [Re £] 2 is shown 
in Fig. 14.21 on the whole energy-momentum plane, and in Fig. 14.31 at fixed mo- 
mentum (in the middle of the grid) of k = 1 GeV as a function of energy. It is 
larger in the chiral limit as in the case with explicitly broken chiral symmetry, 
because [ReE] 2 ~ A 2 ~ (m 2 — m 2 ) 2 is maximal for m 2 = 0. Note that [ReS] 2 is 
small compared to the (squared) 4-momentum independent mass of the u-meson, 
shown in Fig. 14.11 I come back to the influence of [ReS] 2 in the discussion of the 
spectral density. 



w(MeV) 




(4.27) 
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Figure 4.4: The decay width of the a-meson T a , as a function of temperature T and 
a-meson vacuum mass m a , at momentum k = 1 GeV, in the case with explicit chiral 
symmetry breaking (left) and the chiral limit (right). 



The decay width 

After calculating the (4— momentum dependent) imaginary part of the cr-meson 
self-energy (|4.8ajl . the decay width IV (|4.24jl is calculated at the quasiparticle 
energy uj = uo a . The imaginary part of the self-energy, and therefore the decay 
width, of the pion is of order ~ 1/N, and therefore neglected. 
In Fig. 14.41 the decay width of the u-meson IV is shown as a function of temper- 
ature T. The qualitative behaviour is similar in all cases, but the decay width 
is larger in the chiral limit (m n = 0) compared to the case with explicit chiral 
symmetry breaking (m, 7^ 0), because IV ~ ImE ~ A 2 ~ (m 2 — m 2 ) 2 < m^ 
[cf. Eqs. (|4.8aj) and (jl.23|) ]. The decay width is a strictly monotonic decreas- 
ing function with temperature, and becomes approximately zero (equal zero) for 
temperatures larger than ~ 200 MeV (T > T x ) in the case with explicitly bro- 
ken chiral symmetry (the chiral limit). This is a consequence of the (partial) 
restoration of the chiral symmetry, the masses of the chiral partners become (ap- 
proximatively) degenerate, and therefore the phase space of the a — » 2n decay 
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Figure 4.5: The decay width of the cr-meson T a , in the improved Hartree (triangles) and 
in the improved Hartree-Fock (squares) approximation, as a function of the temperature 
T, at momentum k = 1005 MeV and a-meson vacuum mass of m a = 600 MeV, in the 
case with explicit chiral symmetry breaking. 



is squeezed. The dependence on the vacuum mass of the cx-meson is significant. 
The reason for this is that T a ~ ImS ~ m^ which agrees reasonably with the 
results, Im E = 1 : 3 : 16 for m CT = 400 : 600 : 800 MeV. 

In Fig. 14.51 the decay width of the cr-meson in the improved Hartree (triangles) 
is compared with the improved Hartree-Fock |RRR05| (squares) approximation. 
The main difference comes from the combinatorial factor in front of the two- 
pion term (~ / PP) in the imaginary part of the self-energy [cf. Eqs. ()4.6a|) 
with (|4.8a|) ]. A part of this contribution is of order ~ 1/iV, and neglected in 
the Hartree approximation. Therefore, this factor is 2 ■ iV = 8 in the improved 
Hartree, and 2 • (N — 1) = 6 in the improved Hartree-Fock approximation, which 
would lead to a decay width which is an factor pa 1.33 larger in the improved 
Hartree approximation. The remaining difference, shown in the plot, comes from 
the a-meson term (~ 3 ■ 3! / SS) in Eq. f|4.8aj) . which vanishes in the large- N 
limit. 
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Figure 4.6: The spectral density of the a -meson p a , as a function of energy u, momen- 
tum k, and cr-meson vacuum mass m a at a temperature of T = 100 MeV, in the case 
with explicit chiral symmetry breaking (left) and in the chiral limit (right). 



The spectral density 

Finally, after solving the coupled condensate and Dyson- Schwinger equations self- 
consistently, Eqs. (|4.3jl and (j4.5j) . the spectral density of the cr-meson p a is given 
in the chiral broken phase (a ^ 0) by Eq. (J4.15|) . and in the restored phase (a = 0) 
by Eq. (jOBjl . 

In Fig. l4.6l the spectral density is plotted as a function of energy and momentum, 
and in Fig. 14. 71 at fixed momentum (as for the real part in the middle of the grid) 
of k = 1 GeV. In Fig. 14.71 all possible parameter sets (cf. table 4.1) are compared. 
The spectral density does not exhibit a pronounced peak at the mass-shell energy 
uj a = yj k 2 + M%. The reason for this is that the energy of the a-meson u a is large 
enough to decay into two pions. As discussed above, the decay width of this 
process is large and becomes even larger for increasing m CT , as shown in the plot. 
A remarkable difference between the chiral limit (right) and the case with explicit 
chiral symmetry breaking (left) is that p a w around to ~ 1 GeV in the case with 
explicit symmetry breaking but not in the chiral limit. This effect can be traced 
back to the threshold energy for a — > 2tt, which is u = 2M V1 and therefore zero 
in the chiral limit but not in the case with explicit chiral symmetry breaking. For 
uj < k the a-meson is Landau-damped. 
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Figure 4.7: The spectral density of the a-meson p a , as a function of energy to, and o- 
meson vacuum mass m c , at a momentum of k = 1 GeV and a temperature ofT = 100 
MeV, in the case with explicit chiral symmetry breaking (left) and in the chiral limit 
(right). 




Figure 4.8: The spectral density of the a-meson p a , as a function of energy u, at a 
momentum of k = 1 GeV, and a vacuum a-meson mass of m a = 600MeV, in the case 
with explicit chiral symmetry breaking. The results are shown for the improved Hartree 
(full line) and the improved Hartree-Fock (dotted line) approximation at temperatures 
ofT = 160 MeV (left), and 320 MeV (right). 
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To illustrate the effects of neglecting terms of order ~ 1/-/V, I compare in Fig. 14.81 
the results from the improved Hartree (full lines) with the results from the im- 
proved Hartree-Fock |RRR05j (dotted lines) approximation. The results are 
shown in the low-temperature regime, T = 160 MeV (left), and in the high- 
temperature regime, T = 320 MeV (right). In the improved Hartree-Fock ap- 
proximation one does not neglect the ~ 1/N terms in the imaginary part of 
the self-energy, which leads to a nonzero width of the pion spectral density, and 
therefore to a washed-out cx-meson spectral density. 

Another major difference between the improved Hartree and Hartree-Fock ap- 
proximation is that the 4-momentum dependent real parts are neglected in the 
latter one. However, as mentioned, the differences in Fig. 14.81 can be mainly 
traced back to the the neglected terms of order ~ 1/N. To quantify the differ- 
ence arising from the 4-momentum dependent real parts, I calculate the relative 
change between the spectral density with and without [ReS] 2 , averaged over the 
energy-momentum grid, (diff) = (|(pw — Pwo)/pw\)k,w, where p w is the stan- 
dard spectral density given by Eq. (j4.15|) . p wo the spectral density given by 
Eq. (j4.15j) with [ReS] 2 = 0, and (. . .)k,w denotes the average over the energy- 
momentum grid. This mean value is (diff) = 5.61 ± 2.90% for T = 160 MeV, and 
(diff) = 0.44 ± 0.01% for T = 320 MeV. 
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Conclusions & Outlook 



In the following I briefly discuss the major results from the last three chapters 
and explain what is planned in the future. 

Chapter II: The quark mass dependence of the transition temperature 
[DRR04J 

Three-colour QCD exhibits a (weakly) first-order deconfining phase transition 
at a temperature T c j^Ja ~ 0.63 in the limit of infinitely heavy quarks (\/a ~ 
0.425 GeV denotes the string tension at T = in this theory). Near T c , the screen- 
ing mass for the fundamental Polyakov-loop £ drops substantially |KKLL00] , and 
so one might hope to capture the physics of the phase transition using some 
effective Lagrangian for £ |YS82l IFS021 IDFoTl IDP02bl IDP02al ISDJ01I IUM001 
IMMO02I ImSTO! IFuk04j . 

For finite quark masses, a term linear in I appears which breaks the Z(3) center- 
symmetry explicitly. This reduces the deconfinement temperature, with AT C /T* 
on the order of the expectation value of the Polyakov-loop at T~, cf. Eq. (j2.13|) . 
At some point then, the line of first-order deconfinement phase transitions 
ends |GGP94I IGK84I IMU95I lA+99l lB+90| . I have provided a quantitative esti- 
mate of this point, m n ~ 4.2 sja « 1.8 GeV(= m q w 0.9 GeV) and T c ~ 240 MeV 
for Nf = 3 degenerate flavours, by matching the effective Lagrangian for the 
Polyakov-loop to lattice data on T c (m n ) |KLP01j . Assuming that bi oc Nf |A + 99j 
shifts "D" to m n ~ 1.4 GeV for N f = 2 and to 0.8 GeV for N f = 1. 
Going to even smaller quark (or pion) masses leaves a crossover between the 
low-temperature and high-temperature regimes of QCD. The dependence of the 
crossover temperature T c on the pion mass appears to be well described by a 
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small explicit breaking of the Z(3) center symmetry, 61 ~ exp(— m w ), down to 
m Tr/y/< J — 1) which is the smallest pion mass covered by the lattice data of 
Ref. |KLP01j . On the other hand, a linear sigma model leads to a stronger 
dependence of T c on than seen in the data. 

In turn, in the chiral limit, and for Nf = 3 flavours, one expects a first-order 
chiral phase transition |PW84b| IGGP94| IB + 90| . The linear sigma model should 
then be an appropriate effective Lagrangian for low-energy QCD |PW84bt IRVV93 , 
IhRSBOOMMMO P94I lM(7PP92l IGol 831 1( ;( ; P94j . The first-order chiral phase tran- 
sition ends in a critical point "C" if either the mass of the strange quark or those 
of all three quark flavours increase. Given that the explicit symmetry breaking 
term for the Polyakov-loop remains rather small when extrapolated to m n — > 0, 
that is b± — > 0.2, I speculate that "C" might be rather close to the chiral limit. 
Indeed, recent lattice estimates for Nf = 3 place "C" at ~ 290 MeV [KXSOlj 
for standard staggered fermion action and N t = 4 lattices; improved (p4) actions 
predict values as low as ~ 67 MeV jK+04aj . 

Of course, the question arises why, for pion masses down to ~ 400 MeV, the 
QCD crossover is described rather naturally by a slight "perturbation" of the 
= 00 limit, in the form of an explicit breaking of the global Z(N C ) symmetry 
for the Polyakov-loop. Physically, the reason is the flatness of the potential for 
£ in the pure gauge theory at T c , see e.g. the figures in |DP0H IDP02b[ ISDL02j . 
which causes the sharp drop of the screening mass for i near T c + |KKLL00] . This 
is natural, given that finite-temperature expectation values of Polyakov-loops 
at N c = 3 are close to those at N c = 00 |DHL + 04] . where the potential at T c 
becomes entirely flat j( ; W801 IKSS821 IDHL+04] . Hence, a rather small "tilt" of the 
potential (due to explicit symmetry breaking) quickly washes out the deconfining 
phase transition of the pure gauge theory, and causes a significant shift AT C of the 
crossover temperature already for small b\. If so, then for N c — > 00, at the Gross- 
Witten point, the endpoint "D" should be located at 61 = 0; the discontinuity 
for the Polyakov-loop at T c , which in a mean-field model for the pure gauge 
theory is 1/2 at iV c = 00 jDHL+041 IHW8H1 IKSS82] . then vanishes for arbitrarily 
small explicit symmetry breaking. This has previously been noted by Green and 
Karsch |GK84| within a mean-field model. If confirmed by lattice Monte-Carlo 
studies, one might improve the understanding of the degrees of freedom driving 
the QCD crossover for pion masses above the chiral critical point "C". 
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Chapter III: The improved Hartree-Fock approximation [RRR05J 

In chapter III, I have studied the O(N) linear sigma model at nonzero tem- 
perature within a self-consistent many-body resummation scheme. This scheme 
extends the standard Hartree-Fock approximation by including nonzero decay 
widths of the particles. In the standard Hartree-Fock approximation, the self- 
energies of the particles consist of tadpole diagrams which have no imaginary 
part. Consequently, all particles are stable quasi-particles. In order to obtain a 
nonzero decay width, one has to include diagrams in the self-energy, which have 
a nonzero imaginary part corresponding to decay and, in a medium, scattering 
processes. 

In order to incorporate the nonzero decay width in a self-consistent way, I ap- 
ply the Cornwall- Jackiw-Tomboulis formalism. The standard Hartree-Fock ap- 
proximation is obtained by considering only double-bubble diagrams in the 2PI 
effective action, leading to the (energy- and momentum-independent) tadpole 
contributions in the 1PI self-energies. In order to extend the Hartree-Fock ap- 
proximation, I additionally take into account diagrams of sunset topology in the 
2PI effective action. This has the consequence that the 1PI self-energies obtain 
additional energy- and momentum-dependent one-loop contributions which have 
a nonzero imaginary part. The spectral densities of a- mesons and pions are then 
computed as solutions of a self-consistent set of Dyson-Schwinger equations for 
the cr-meson and pion two-point functions, coupled to a fix-point equation for 
the chiral condensate. I only take into account the imaginary parts of the new 
one-loop contributions. I made sure that the spectral densities obey a standard 
sum rule by adjusting their normalisation, if necessary. 

I found that the temperature T x for chiral symmetry restoration is about 20% 
smaller as compared to the Hartree-Fock approximation when including nonzero 
particle decay widths. My value T x ~ 175 MeV agrees reasonably well with 
lattice results |LP03| . I computed the decay widths of cr-mesons and pions as a 
function of temperature. The vacuum value for the cr-meson decay width comes 
out to be in the experimentally observed range, without adjusting a parameter 
of the model. It stays approximately constant up to temperatures ~ T x and then 
decreases sharply with temperature. The pion decay width grows from zero at 
T = to a value ~ 100 MeV at T ~ T x , and then also decreases with temperature. 
I also investigated the spectral densities of cr-mesons and pions as functions of 
energy and momentum for temperatures in the range of T = 80 to 320 MeV. 
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Below the chiral phase transition, the spectral density of the cr-meson is broad, 
due to the possible decay into two pions. It develops a peak at the quasi-particle 
mass shell above the chiral phase transition, when this decay channel is closed. 
On the other hand, the spectral density of the pion always exhibits a distinct peak 
at the quasi-particle mass shell. The width of this peak is due to scattering off 
cr-mesons in the medium. Above the chiral phase transition, the spectral densities 
of cr-mesons and pions become degenerate in shape. 

Chapter IV: The improved Hartree approximation |R6d05j 

In chapter IV, I systematically improved the standard Hartree (or large- N) ap- 
proximation of the O(N) linear sigma model by taking into account, additionally 
to the double-bubble diagrams, the sunset diagrams in the effective potential 
of the CJT formalism. This leads to 4-momentum dependent real and imagi- 
nary parts of the Dyson-Schwinger equations. In contrast to chapter III, I didn't 
neglect the 4-momentum dependent real parts, to study their influences to the 
results. I solve these and the equation of the condensate selfconsistently in the 
chiral limit (m u = 0), and in the case with explicitly broken chiral symmetry 
(m n 7^ 0), to get the decay width and the spectral density of the cr-meson. 
First, I presented the results for the real parts of the Dyson-Schwinger equations 
and the condensate equation. The 4-momentum independent parts, i.e., the ef- 
fective masses and the chiral condensate exhibit a crossover transition in the case 
with explicit chiral symmetry breaking and a second-order phase transition in 
the chiral limit, which agrees with the predictions made by Pisarski and Wilczek 
PW84aJ. I found, that the 4-momentum dependent real part of the self-energy 
is rather small compared to the (squared) 4-momentum independent effective 
masses. It is larger in the chiral limit, simply because [Re] 2 ~ A 2 ~ (m 2 — m 2 ) 2 
is maximal for m 2 = 0. 

The decay width shows qualitatively the same behaviour in all cases. It is a 
decreasing function with temperature and becomes (approximatively) zero in the 
high-temperature regime. Nevertheless, quantitatively it depends strongly on the 
choice of the vacuum mass of the cr-meson, by reason that r a ~ Im X ~ A 2 ~ m*. 
Also, the decay widths are larger in the chiral limit, due to the fact that T a ~ 
A 2 ~ (m 2 - ml) 2 . 

In the low-temperature regime, the spectral density of the cr-meson is a very broad 
function in energy (and becomes even broader for larger vacuum mass m a ) due to 
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the a — > 27r decay. In the high-temperature regime, the spectral density becomes 
more and more a delta function. Remarkable is the fact that p a becomes zero in a 
certain energy interval, in the case with explicit chiral symmetry breaking but not 
in the chiral limit. This can be traced back to the threshold energy for o — > 2tt, 
which is 2m n , and therefore zero in the chiral limit but not in the other case. 
For uj < k the u-meson is Landau-damped. The influence of the 4-momentum 
dependent real part of the self-energy to the spectral density is rather small and 
does not change the results qualitatively. 
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Outlook 

The present studies can be continued along several lines: 

• Throughout this thesis I used the linear a model with 0(4) symmetry, which 
is isomorphic to the linear a model with SU(2) x SU(2) symmetry. Both 
models contains only mesons composed from the lightest two quarks (the up 
and down quarks), but the 0(4) model is the limit of the SU(2) x SU(2) 
model with maximally broken U(1)a symmetry, which leads to the fact 
that the r\ and the mesons become infinitely heavy, and hence cannot be 
describted by this theory, cf. RRR03J. It would be exciting to study the 
influence of nonzero decay width effects also for these mesons, and addition- 
ally for mesons with strange and maybe also charm degrees of freedom. The 
scalar (s) and pseudoscalar (ps) mesons contained in the linear a models 
are summarized in the following table: 





0(4) 


SU{2) r x SU(2) e 


SU(3) r x SU(3) t 


SU(4) r x SU{4) e 


s 


a 


n ± n° 


kP 2 , kP, K°, /o 


Do, Dq,Dq, Df , Xc,o 


ps 




V 


K ± ,K°,K°,ri' 


D ± , D°, D°, Df, t) c 



In chapter II, the chiral critical point is at m q = 0, because only the 2-flavor 
case is discussed. Extending the model by the strange degree of freedom 
would allow us to give an estimate for the position of this point in the whole 
quark mass plane, discussed in Sec. 1.3, cf. Fig. 11.21 

So far only the scalar and pseudoscalar mesons are taking into account. The 
inclusion of baryonic degrees of freedom into a framework, similar to that 
discussed in Chapters III or IV, is still under investigation |Bec05j . Also, 
the inclusion of vector mesons [RRQU IStej is under investigation, and is of 
particular importance, since in-medium changes in the spectral properties of 
vector mesons are reflected in the dilepton spectrum RR05 which, in turn, 
is experimentally observable in heavy-ion collisions at GSI-SIS, CERN-SPS 
and BNL-RHIC energies. 

As disussed, a major advantage of effective models over lattice QCD cal- 
culations is that no problems with a nonzero chemical potential occur 
|VV()02[ ICOVV05] . Including the chemical potential in the framework 
of the improved Hartree or Hartree-Fock approximation would us allow to 



study, e.g., the influence of nonzero decay width effects on the chiral critical 
endpoint in the T — fi— plane shown in figure PHIl (cf. e.g. SMMROlJ). 



Up to now, I neglected the ultraviolet divergent vacuum parts of the in- 
tegrals, i.e., I just used trivial renormalization. One should check what 
happens if one use a "real" renormalization scheme (e.g. the cut-off or 
counter-term scheme, cf. jLROOj ). 

The influence of the 4-momentum dependent real parts of the self-energy in 
the improved Hartree approximation is very small, nevertheless, one should 
check its influence also in the full Hartree-Fock approximation. 
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-Appendix A- 

The calculation of the 
diagrams 

A.l In the improved Hartree-Fock approxima- 
tion 

In this Appendix I discuss the calculation of the three types of diagrams (see 
Fig. IA.1|) contributing to the self-energies and the condensate equation. To eval- 
uate them explicitly, I use standard techniques of thermal field theory, see for 
example Refs. |l)as971 |Kap93j IT/BOOl lhar.971 lhvW87l Ml ■ 

The tadpole diagram 

In the following I discuss the calculation of the tadpole diagram in Fig. IA.1I a as 
a functional of the spectral density. In the imaginary-time formalism it is given 




Figure A.l: The general topology of the tadpole diagram a, the cut sunset diagram b, 
and the sunset diagram c. 
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The calculation of the diagrams 



by 



(a.i; 



where T is the temperature, and A is either the cr-meson or pion propagator. 
The first step is to perform the sum over the Matsubara frequencies. To this aim 
I introduce the mixed representation of the propagator 



,1/T 

A(-iu n ,q)= dTexp(-iu n T)A(T, q) 
Jo 



with 



/OO fajj 
— [6(r) + f(u)] p(u, q) exp(-^r) , 
-oo AH 



(A.2) 



(A.3) 



where f(u) = [exp(u;/T) — 1] 1 is the Bose-Einstein distribution function. With 
the identity 

oo , 

T^exp(-z^ n r)= 8 (r - 



m 
T 



(A.4) 



the Matsubara sum can be performed analytically 



/d 3 q f^/T 
— — ^ rfrexp(-^ n r)A(r,q) 



f d°q f 1 / 1 , . . 
d 3 q 1 



5{t) + 8[t- 



T 



(2tt) 3 2 
d 3 q 



A(0,q) + A(-,q 



(2tt)* 



A(0,q) 



(A.5) 



where use has been made of the KMS condition A(r, q) = A(r — 1/T, q). With 
Eq. (TOl) I finally have 



r 



00 duo f d 3 q 



2tt J (2ttY 



p(u,q) 



(A.6) 



Due to isotropy of space, the spectral density of a scalar particle cannot depend 
on the direction of q, thus the angular integration can be carried out: 



T = 4tt 



did 



2tt Jo (2tt) 



dq 2 
— 9 



+ /(") 



p(u,q) 



(A.7) 
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Using the fact that the spectral density for bosonic degrees of freedom is an odd 
function of the energy, p(u) = —p(—u), 

T - *f £^i/m-/<-"M-<.»> 

Subtracting the divergent vacuum contribution I obtain: 

r= (2^/o ^gVM/K^g) ■ (A.9) 
The cut sunset diagram 

In this section I calculate the imaginary part of the cut sunset diagram of Fig. IA.1I 
b. In the imaginary-time formalism the cut sunset diagram is given by 

C(-iw m ,k) = / A 1 (K-Q)A 2 (Q) 

JQ 

= T E / 7 |^A 1 H(o; m -a; n ),k-q)A 2 Ha; n ,q) ,(A.10) 

where A 1; A 2 are the propagators of cr-mesons and/or pions. The major difference 
between the tadpole diagram, discussed in the last section, and this diagram is the 
fact that it explicitly depends on the external four-momentum, = (— iu m , k). 
Analogously to the last section one introduces the mixed representation [see 
Eq. (jA.2j) ] of the propagators Ai and A 2 

— — - / dnexp[-z(a; m -a; n )ri]Ai(Ti,k-q) 
(27r) d Jo 

,1/T 

x / rfr 2 exp(-zc t ; n r2)A2(r2,q) . (A.ll) 
Jo 

The summation over the Matsubara frequencies leads to a delta function [cf. 
Eq. (H3J)] 

/(pq r^/T 
^— ^ J Q dTexp(-iu m T)Ai(r, k- q)A 2 (r, q) . (A.12) 

Introduction of the spectral densities [see Eq. (|A.3|) ] for both propagators in the 
mixed representation leads to 

d 3 q fV T 



C(-iu m ,k) = J j^y J o drexp(-iuj m T 
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x / ^-[ 1 + /(^i)]Pi(^i ) k-q)exp(-^ir) 

J-oo 2lX 

/OO (l(jQn 
— [1 + f(u 2 )}p 2 (u 2 , q) exp(-^ 2 r) . (A.13) 
-oo Z7T 

The integration over r can be performed analytically. Employing the identity 

exp - l] [1 + [1 + /(u*)] = -1 - fM - f(u 2 ) , (A.14) 

the expression can be rewritten as 
n( ■ i \ f°° dwi du 2 1 + fM + fM f d 3 q 

C{-iu m , k) = / — : ■ ■ / — — Pl {ui, k - q)p 2 (w 2 , q) • 

J-oo 2ir 2tt tuj m + uji + uj 2 J {2it) 6 

(A.15) 

Using the fact that the spectral densities do not depend on the direction of 
momentum and 

X : f(\y - -x\)g(x) = —jj-^f dx 1 x 1 dx 2 x 2 f(xi)g(x 2 ) 



(2tt) 3JXW j y (2vr) 3 7o 

x (\xi - x 2 \ < y < xi + x 2 ) , (A.16) 

one can perform one angular integration, yielding 

dujt du 2 1 + f(u x ) + f{oj 2 ) 



C(-iuj rn , k) = J 



oo 27T 2lT iu m + IjJ\+ u 2 

2tt 1 r°° 
x -JTJ^js J Q d( li ( li d( l2q2pi(u 1 ,q 1 )p 2 (uj 2 ,q 2 ) 

X e{\qi-q 2 \ <k< qi +q 2 ) . (A.17) 

The imaginary part of the retarded cut sunset diagram can be extracted using 
analytical continuation, —iu) m —>■ u + is, and the Dirac identity Iml/(x + ie) = 

-7r5(x), 

lmC(u,k) = ^^-j_^du)xdu 2 [l + f{u x )+f{u 2 )]5{u-ux-u 2 ) (A.18) 



X 



2(2vr) 

POO 

/ dq 1 q i dq 2 q 2 Q (\q x - q 2 \ < k < q x + q 2 ) Pi(u u q x ) p 2 (co 2 , q 2 ). 
Jo 

This integral is finite and therefore does not require renormalization. 
The sunset diagram 

In the imaginary-time formalism the sunset diagram shown in Fig. IA.1I c is given 
by the following expression: 

/d^l d?q 
(27r)3 (2tt}3 Al (~^ m > l) A 2(-^ n , q)A 3 (-z(-^ n - u m ), -q - 1) . 

(A.19) 



run 
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I introduce the mixed representation, see Eq. (|A.2jl . for the three propagators 
Ai, A 2 , and A 3 . One can perform the Matsubara sums employing Eq. (jA.4j) . In 
the mixed representation one introduces the spectral densities pi, p 2 , and p 3 in 
order to perform the r integration analytically. Employing the identity 

exp (_ ^+^ + ^ - l] [1 -h fMl [1 + fM] [1 + fMl 
= -1 - fM - fM - fM - fMfM - fMfM - f(ou 2 )f(co 3 ),(A.20) 

for the Bose-Einstein distribution functions I obtain 

duii du 2 du 3 
—oo 2n 2n 2tt 

1 + fM + fM + f(u 3 ) + fMfM + fMfM + fMfM 

X 

LO\ + UJ 2 + ^3 

[ d 3 l d 3 q 

x J 7^p^pPi(^i,l)P2(^2,q)p3(^3,-q-l) ■ (A.21) 

Using the fact that the spectral densities do not depend on the direction of 
momentum, the angular integrals can be carried out with the result 

2 roo roo 

S = 7n~Y7 / dcu 1 duj 2 duj 3 dq l q 1 dq 2 q 2 dq 3 q 3 Q(\q 1 -q 2 \<q 3 <q 1 + q 2 ) 

(2lT)< J-oo J 

v i + fM + fM + fM + fMfM + fMfM + fMfM 

x 

+ CJ 2 + ^3 

x PiMQiiPiMOadPaMQa) ■ (A.22) 
Using the antisymmetry of the spectral densities, I obtain 

4 roc 

S = j—r^ du 1 duj 2 duj 3 dq 1 q 1 dq 2 q 2 dq 3 q 3 6 (\qt - q 2 \ < q 3 < q 1 + q 2 ) 
(2n)' Jo 

fMfM + fMfM + fMfM + fM + fM + fM + 1 



X 



L} 1 + 0J 2 + <JJ 3 

| fMfM - fMfM + fMfM + fM 

U)\ + UJ 2 — UJ 3 

| fMfM - fMfM + fMfM + /M 

CJi — W 2 + UJ 3 

+ fMfM + fMfM - fMfM + fM 



—U)\ + UJ 2 + UJ 3 

x Pi(wi,gi)p 2 (w2, 92)^3(^3, 93) • (A.23) 



<s<s 



The calculation of the diagrams 



Finally, I discard the ultraviolet-divergent parts, and relabel the integration vari- 
ables 

4 roc 

$ = 7n~\7 i du 1 du 2 duj 3 dq 1 q 1 dq 2 q 2 dq 3 q 3 Q (\qx - q 2 \ < q 3 < Qi + fe) 
{2tt)' Jo 

X fMf(u 2 )( ■ 1 ■ + J— + — + 1 



UJi + Co> 2 + Co> 3 — CJi +0^2 + ^3 UJi — UJ 2 + UJ 3 —UJ\ — UJ 2 + UJ 3 
X [pi(u;i, gi)p 2 (^ 2 , <?2)P3(^3, Qz) + Pl(^l, gi)p 2 (^3, <?2)P3(^2, qz) 

+pi(uj 2 ,q 1 )p 2 (uj 3 ,q 2 )p 3 (uj 1 ,q 3 )] . (A.24) 



A. 2 In the improved Hartree approximation 

In contrast to the improved Hartree-Fock approximation, the spectral density 
of the pion, p n , in the improved Hartree approximation is just a delta-function, 
because the imaginary part of the pion self-energy if of order ~ 1/N and thus 
neglected. Therefore the equations for the tadpole and the cut sunset diagrams, 
Eqs. (jA.9|) and (|A.17|) simplify, but one has to calculate additionally the real part 
of the cut sunset diagram. The full sunset diagram, discussed in IA.1[ is of order 
~ 1/N and is therefore neglected in the large- N limit. 



The tadpole diagram 

The Dyson-Schwinger and the condensate equations of the improved Hartree 
approximation contain exclusively the spectral density of the pion. Therefore one 
has to solve equations of the form 

4 poo 

with different combinatorial factors in front. (Note that this is the version where 
the divergent term is already neglected, i.e., the equation in the trivial renormali- 
sation scheme.) Using the fact that the spectral density of the pion is just a delta 
function, p n (uj, k) = ir/ui w (k){5[ui — u n (k)] — 5[uj + u> n (k)]}, one gets 

T= 2^i ^<7 2 K(<7)]-7K(<7)]. (A.26) 
The cut sunset diagram 

In this section I derive the equations for the 4-momentum dependent imaginary 
and real part of the self-energy (arising from the sunset diagram). As discussed 
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in Sec. A.l this diagram can be expressed as a function of the spectral density 
p n in the imaginary time formalism 

E{-iu m ,k) = ttttit / dcJidcj 2 — ■ ■ A.27 

\ N J (27T) 4 k J-oo lUJ m + UJ\ + u> 2 

POO 

x / dqiqx dq 2 q 2 O (\qi - q 2 \ < k < q x + q 2 ) p^i, qx)p 7r (u 2 , q 2 ), 
Jo 

where f(u>) = l/[exp(o;/T) — 1] is the Bose-Einstein distribution function, and 
u m = 2nmT are the Matsubara frequencies. To extract the imaginary part of 
the retarded self-energy, one uses the Dirac identity, Iml/(x + ie) = — ir5(x), and 
analytic continuation, —ioo m — > uj + ie, 

(A\a\ 2 N 1 r°° 
ImS(w,k) = \-jfJ 2(27r) 3 k J-oo ^ dUJ2 ^ + ^ Ul ^ + ^ U2 ^ ^ 

roo 

x / dq x q l dq 2 q 2 9 (|g x - q 2 \ < k < qi + q 2 ) p- K {uJi,q l )p- K {uo 2 , g 2 ).(A28) 

JO 

As discussed in the improved Hartree approximation of the O(N) model the 



spectral density of the pion is just a delta function (j4.12j) . 
f AXa\ 2 N 1 r°° 

ImE(w,k) = — — —7/ du 1 du2[l + fM + f(u} 2 )]6(u -U1-U2) 

\ iV y 07T fc J-00 
/•oo 

x / dq 1 q 1 dq 2 q 2 9(|gi - q 2 \ < k < q x + q 2 ) [^(gi)^^)] -1 
Jo 

x {5[ui - uj^qi)) 5[lj 2 - u ir (q 2 )] + 5[ui + uv(gi)] 6[u 2 + uj n (q 2 )] 
- 8[ui - Unfa)] 8[u 2 + uj n (q 2 )] - 8[ui + ^(gi)] 8[u 2 - v n (q 2 )]} , 

(A.29) 



where tcv(g) = yg 2 + M 2 is the quasiparticle energy of the pion. To simplify this 
expression one uses the 8{u — uji — oj 2 ) function to carry out the ^-integration, 

/4Acr\ 2 N I r 00 

roo 

x / dqiqidq 2 q 2 Q{\qi-q 2 \<k<qi + q 2 )[u 1T (qi)uj v {q2)Y 1 
Jo 

x {S[ui - u w (qi)] 8[uj -Ui- u w (q 2 )] 
+8[ut + 5 [a; - Wj + 0^(92)] 

-<5[^i - ^(gi)] 8[uj + uj n (q 2 )} 

-Slux + ^(gi)] 6[u - ui - ^(g 2 )]} , (A.30) 
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and the Wi-integration with the help of the 5[uji ± ^(gi)] -functions 



Im Tj(uj, k) 



(A.31) 



= [if) ^kJo ****** 

x © (|?i - 92 1 < k < qi + q 2 ) [u^gx)^^)] -1 
x ( I 1 + /Kfai)] + /[<*> - u v (qi)]}6[w - uj w ( qi ) - u^(q 2 )) 
+{1 + f[-ouM] + f[u + uM]} 5[uj + uM + uo^q 2 )} 
-{1 + /[^(gi)] + /[w - ^(gx)]} <5[u; - ^(gx) + cj 7r (g 2 )] 
-{1 + /[-^(gi)] + f[u + ^(gi)]} 5[o> + ^(gi) - cj 7T (q 2 )]) . 



The 5[u + u n (qi) + u n (q 2 )] function has no support, because ui > and u n (q) > 0. 
Note that evaluating the 5[u — u n (qi) + u n (q 2 )} and the 5[u + u) n (qi) — dj 7r (g 2 )] 
function, the latter two terms cancel each other. The remaining delta function 
can be transformed to a delta function in momentum-space, 



^(go) 



[%i - 9o) + %i + go)] 



where 



q ^J[u-u^q 2 )f-Ml 



(A.32) 



(A.33) 



is the root of the argument of the delta function. Note that in Eq. (jA.31|) q\ > 
and therefore there is no support of the S(qi + go) function. Carrying out the 
g!-integration, using u — u n (q ) = uj n (qi), and relabelling q = q 2 , one gets 



Im £(co>, k) 



'4Aa\ N 1 r 00 



\ N J 8irkJo 
x {1 + /K(g )] + /K(g)]}. 



POO 

/ dq g[w 7r (g)] _1 ©(|g -g| < k < g Q + g) 
Jo 

(A.34) 



To calculate the limit k — * 0, one starts best with Eq. (jA.31|) . uses the following 
transformation 



l + /K(9i)] + /[^-^(gi)] 



1 — exp ( — 



w 
T 



(A.35) 



x {l + /K(g 1 )]}{l + /[o;-a; 7r (g 1 )]} 



the identity 



@(|gi ~qi\ <k<qi + q 2 ) 
hm = 2 8( gi - q 2 ) 



(A.36) 
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(to perform the gi-integration), and relabels q = q 2 , to get 



(AXa\ N 



2 

1 - exp (-- 



dqq 2 5[u-2^{q)][ul{q)}- 1 



x {l+/M«)]}{l + /[w-u,(g)]}. (A.37) 

As explained in |Ris98j this expression can be calculated analytically, using the 
Spence integral, 



Note that in |Ris98j the calculation is performed for the special case uj = m CT , but 
this can be generalised without further problems. 

To calculate the real part of Eq. (jA.27|) . one has to integrate over the principal 
value of + u 2 — u), denoted by P / . . ., 

ReE a;,k 2 = — - 77^t p / dwi<Jw 2 A.39 

X / dq 1 q 1 dq 2 q 2 © (|<?i -^l < fc < 9i + fe) Ptt(^i, 9i)Ptt(^2, 92)- 

JO 

Again, one uses the fact that the spectral density of the pion is just a delta 
function to perform the u\- and the ^-integration, and employs trivial renormal- 
isation, i.e., neglects the divergent parts. This leads to 

/4AoV N 1 r°° 
[Re£(w,k)] 2 = (-— j ^PJ Q dq iqi dq 2 q 2 (A.40) 

x © (|?i - 92 1 < < 9i + 92) [u} n (qi)tj n (q 2 )]~ 1 
x f /K(gi)] + f[w v {<ti)] /K(9i)] + /Kfe)] 



/K(9i)] + , -/K(9i)] + /Kfe)] 



Wtt(9i) - ^(92) - u LU^qt) - uj n (q 2 ) + lu J 

To evaluate the principal value numerically, in an appropriate way, one per- 
forms the following steps. First one introduces a new variable x = u} n (qi) and 
transforms the gx-integration to an x-integration, with dq\ = uj^qij/qidx, and 



qi = \Jx 2 — M 2 = q*. Relabelling q = q 2 leads to 

[ReE(o;,k)] 2 = ~P JJ' dgq^q^dxe (\^ -q\<k<^ + q) 
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'/(*)+ /KG?)] , /(*) + /K(ff)] 



+ 



x + u; 7r (g)— u; x + uv(g) + u; 
-/W + /W9)] -/(*)+/[«*(?)] 



x — u>Tr(q) — id x — ui n (q) + 10 



(A.41) 



Second one transforms the x-integration to a sum over Xi, and uses the mean- 
value theorem to put the terms with the distribution functions in front of the 
integrals 

[Re£(u;,k)] 2 = f-^J dqqMq)]- 1 ^ (\q* - q\ < k < q* + q) 



x + u w (q) - io 



+{/(*) + 

+{-/(*) + /Mff)]} 



+ {-f(Ss) + f[ Uw (q)}} 



dx 



x + uj^(q) + lu 



dx 



x — LOn(q) — to 



dx 



x - uj w (q) + uj 



(A.42) 



where Xi < x < x i+ i, and q* = Jx — M|. Finally, the x-integrals can be per- 
formed analytically, 



[Re£(w,k)] 2 



'4AoV N 1 r°° 



\ N J 8n 2 k 
x ({/(x)+/K(g)]}ln 



poo 

/ dqq^q)}- 1 ^® {\q* ~ q\ < k < q* + q) 
Jo j 

K(g) + a^+i] 2 - ^ 2 



+ {-/(x) + /K(g)]}ln 



w v (q) + Xi] 2 - uj 2 



In the limit k — > 0, using Eq. ()A.36|) . the g-integration can be carried out, 

[ReE( W ,o)] a = (jj^E^wor 1 

K(g*) +x i+1 } 2 -uj 2 



(A.43) 



(A.44) 



x {/(x) + /K(g*)]}ln 



{-/(*) + /K(g*)]} In 



- x i+ if - uj 2 



[w„(q*) - x^ 2 - uj 2 
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